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1. Introduction 

High-speed laminar-to-turbulent transition and turbulence affect the control of flight 
vehicles, the heat transfer rate to a flight vehicle’s surface, the material selected to protect such 
vehicles from high heating loads, the 
ultimate weight of a flight vehicle due to 
the presence of thermal protection systems, 
the efficiency of fuel-air mixing processes 
in high-speed combustion applications, 
etc. Gaining a fundamental understanding 
the physical mechanisms involved in the 
transition process will lead to the 
development of predictive capabilities that 
can identify transition location and its 
impact on parameters like surface heating. 

Currently, there is no general theory that 
can completely describe the transition-to- 
turbulence process. However, transition 
research has led to the identification of the 
predominant pathways by which this 
process occurs. Figure 1.1, taken from Ref. 

1, graphically depicts the known pathways 
by which by which an initially stable, 
laminar boundary layer can become 
unstable, undergo a process of instability 
growth, and eventually break down into 
turbulence. For a truly physics-based model of transition to be developed, the individual stages in 
the paths leading to the onset of fully turbulent flow must be well understood. This requires that 
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Figure 1.1: Transition pathways from laminar to 
turbulent flow. Image taken from Ref. 1. (Reprinted 
with permission of the author and the American 
Institute of Aeronautics and Astronautics). 
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each pathway be computationally modeled and experimentally characterized and validated. This 
may also lead to the discovery of new physical pathways. 

The path to transition in Figure 1 begins with the presence of an initial disturbance, either in 
the freestream flow or within the boundary layer itself. Disturbances can be classified into one of 
three types: entropy spots (related to temperature fluctuations), vorticity, and sound. 
Disturbances in the freestream can originate in the settling chamber of a high-speed wind tunnel 
and pass into the test section. Entropy spots and vorticity originating in the settling chamber can 
be convected into the test section. 3 ’ 4 Sound disturbances can originate in the settling chamber, 
from the wind tunnel nozzle walls, and along the wind tunnel walls. 3,4 Sound radiating from the 
nozzle walls is of particular concern as it is the dominant source of noise in high-speed wind 
tunnels, 5 with its magnitude scaling with Mach number to the fourth power. 6 These disturbances 
can significantly alter the transition-to-turbulence process. 6 ' 9 Low-disturbance, high-speed wind 
tunnels have been developed to mitigate these disturbances. 5 These quiet tunnels provide a 
hypersonic freestream flow with disturbance levels similar to those encountered in flight, with 
freestream pressure fluctuations of less than 0.1%. 6 However, they are limited in extent of 
Reynolds number, Mach number, and enthalpy operability. 10 Table 1.1 provides velocity (U), 
temperature (T), mass-flow (m), and pressure (P) fluctuation levels for two conventional 
hypersonic tunnels 10 ' 12 and one quiet hypersonic tunnel 13,14 as measured with hot-wire (U, T, m ) 
or Pitot (P) probes. 


Table 1.1: Measured or computed maximum fluctuation magnitudes in high-speed wind 

tunnel facilities. 


Tunnel 

Freestream 

Property 

Von Karman Institute Mach 6 
Hypersonic Wind Tunnel H3 

NASA Langley 
20-lnch Mach 6 

Purdue Boeing/ AFOSR 
Mach 6 Quiet Tunnel 

U’/U 

0.6 %* (Ref. 15) 

- 

- 

T'/T 

0.8 % (Ref. 15) 

0.20 % (Ref. 16) 

- 

m'/m 

5.2 % (Ref. 15) 

1.01 % (Ref. 16) 

0.2 %*" (Ref. 13) 

P'/P 

1.0 % (Ref. 15) 

1.6 % (Ref. 17) 

0.01 % (Ref. 18) 


Derived value 

Measurement made within laminar boundary layer 


Disturbances that originate within the boundary layer can be caused by a number of factors, 
including the presence of discrete two-dimensional and three-dimensional roughness elements, 
cavities, distributed roughness, and mass injection. Various roughness types, cavity types, and 
mass injection schemes have been used in transition experiments at NASA Langley. ' The 
influence of roughness on the hypersonic transition process varies, although transition typically 
moves forward (i.e. upstream) with increasing roughness height. A roughness Reynolds 
number of the fonn Re k — p k U k k/ ji k , with k being the roughness height and all other 
parameters (density, p, velocity, U, and viscosity, p) evaluated at k, is often used in correlation- 
based approaches to estimate transition location or the effectiveness of inducing transition 
behavior in the wake region or immediate vicinity of the trip. There is no physical mechanism- 
based theory for roughness-induced transition, although roughness can generate streamwise 


vorticity, a wake shear layer, and can interact with incident freestream disturbances. These 
disturbances then enter the boundary via the receptivity process Streamwise vorticity generation 
can result in transient growth instability amplification. ’ Mass-flow fluctuation amplitude and 
frequency measurements in the wake of a roughness element in a quiet tunnel were performed in 
Ref. 13, with sample results are listed in Table 1.2. Note that the percentage fluctuations are ~2 
orders of magnitude larger than the freestream fluctuations in a quiet wind tunnel (Table 1.1). 


Table 1.2: Measured or computed maximum fluctuation magnitudes in high 
speed transition and turbulent flows. 


Property 

Transition downstream of a large 
roughness element in a laminar 
boundary layer 

Fully Developed Hypersonic 
Turbulent Boundary Layer 
(non-reacting) 

U'/U 

25 % (Ref. 26) 

15 % (Ref. 27) 

T'/T 

33 % (Ref. 26) 

47 % (Ref. 27) 

m' /m 

43 % (Ref. 13) 

52 % (Ref. 27) 

P’/P 

5 % (Ref. 13) 

7.5 % (Ref. 27) 


Receptivity is the process by which either the freestream or roughness-generated 
disturbances mentioned are coupled into the boundary layer flow, generating and/or enhancing 
instabilities." Returning to Fig. 1, the magnitude of the imparted disturbance determines which 
pathway to turbulence they will follow. Reviews of the various instability types encountered 
along paths A-E in Fig. 1 can be found in the literature. 29 ' 34 Both computational and 
experimental studies have been performed to study the interaction of the disturbances with the 
various instability mechanisms. The majority of experiments used to conduct these studies 
involve physical probes, which are either surface mounted or placed in the boundary layer flow 
itself, to measure disturbance parameters. In low-speed environments, the magnitude of these 
disturbances varies from a few hundredths of a percent to -10%. Similar magnitudes can be 
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expected for higher-speed flows. Experiments characterizing normalized amplitude growth of 
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instabilities for supersonic and hypersonic cases have also been conducted. 

Fully developed turbulent flowfields are computed using several different computational 
fluid dynamics (CFD) approaches including Reynolds Averaged Navier Stokes (RANS), Large 
Eddy Simulation (LES) and Direct Numerical Simulation (DNS) and other methods, including 
hybrid techniques. Unfortunately, experimental high-speed hypersonic turbulence data is very 
sparse. For experimental data in this regime that does exist, much of it is absent chemistry 
effects, which play a critical role in high hypersonic Mach number flows. 40 Even when limiting 
discussion to perfect gas experimentation, measurements of turbulent fluctuations in the 
hypersonic regime with hot-wire probes appear to suffer from many of the limitations, as 
discussed below. 40 Therefore, molecular-based optical diagnostics seem well-poised for both 
non-reacting and reacting high-speed turbulent flow measurements. A survey was recently 
commissioned by NASA to look at simulation requirements, technology gaps and critical 
impediments that need to be overcome to greatly enhance CFD analysis and design capabilities 
by the year 2030. 41 The report found that “The single greatest impediment mentioned by a 
majority of the survey participants is the lack of foundational building-block validation datasets 


for physical model development.” Such datasets would benefit from optical diagnostics such as 
those described herein. The right hand column of Table 1.2 shows computed property 
fluctuations in a fully developed non-reacting boundary layer. 

Tables 1.1 and 1.2 provide indications of the measurement precisions that would be required 
to resolve the relevant flow property fluctuations in different applications. These data provide a 
guide for comparison with the individual measurement technique’s capabilities detailed below. 
Measurement precision is a key consideration in developing measurement technology for 
transitional and turbulent flows because precise instruments are needed to measure the 
fluctuations and the growth of the fluctuations. Ideally, an instrument’s precision would be 
negligible compared to the fluctuations in the parameter being measured. When choosing to 
implement a particular measurement approach, additional considerations are also important. 
These include considerations of the required accuracy, frequency response, the prior experience 
of the research team, available equipment, available budget and other factors. These factors are 
discussed in more detail in the Conclusions. 

As stated previously, some of the more common measurement devices for flowfield 
transition-to-turbulence and turbulent boundary layer measurements include hot-wires and Pitot 
probes. While hot-wire probes can have a response time and dynamic range appropriate for high- 
speed measurements, they only provide data at a point, are difficult to calibrate, are fragile, 
have limited bandwidth, 37 and are intrusive. 10 Pitot probes, while more robust than hot-wires, 
have similar issues. Since the fluctuation levels associated with the onset of instabilities is an 
order-of-magnitude less than those associated with turbulence, the intrusive nature of such 
probes can have an adverse effect on the parameter being measured. Limitations on bandwidth 
preclude these physical probes from resolving high-speed, high-frequency instability modes, 
with bandwidths of ~1 MHz being desired. Additionally, the point-wise nature of data 
acquisition prevents the detailed study of spatial behavior of transition-to-turbulence and fully- 
developed turbulent behavior. A review of some of the available non-intrusive optical diagnostic 
capabilities with applicability to transition-to-turbulence measurements is available in Ref. 37. 

This document is intended to describe molecular based measurement techniques that have 
been developed, addressing the needs of the high-speed transition-to-turbulence and high-speed 
turbulence research fields. In particular, we focus on techniques that have either been used to 
study high speed transition and turbulence or techniques that show promise for studying these 
flows. This review is not exhaustive. In addition to the probe-based techniques described in the 
previous paragraph, several other classes of measurement techniques that are, or could be, used 
to study high speed transition and turbulence are excluded from this manuscript. For example, 
surface measurement techniques such as pressure and temperature paint, phosphor 
thermography, skin friction measurements and photogrammetry (for model attitude and 
deformation measurement) are excluded to limit the scope of this report. Other physical probes 
such as heat flux gauges, total temperature probes are also excluded. We further exclude 


measurement techniques that require particle seeding though particle based methods may still be 
useful in many high speed flow applications. 

This manuscript details some of the more widely used molecular-based measurement 
techniques for studying transition and turbulence: laser-induced fluorescence (LIF), Rayleigh 
and Raman Scattering and coherent anti-Stokes Raman scattering (CARS). These techniques are 
emphasized, in part, because of the prior experience of the authors. Additional molecular based 
techniques are described, albeit in less detail. Where possible, an effort is made to compare the 
relative advantages and disadvantages of the various measurement techniques, although these 
comparisons can be subjective views of the authors. Finally, the manuscript concludes by 
evaluating the different measurement techniques in view of the precision requirements described 
in this chapter. Additional requirements and considerations are discussed to assist with choosing 
an optical measurement technique for a given application. 

2. Laser-Induced Fluorescence (LIF) 

Laser-induced fluorescence (LIF) uses a laser to probe individual species within the 
flowfield, providing information pertaining to both the thermodynamic (pressure, temperature, 
mole fraction) and fluid dynamic state (velocity) of the gas. The laser can either be focused to a 
point for LIF measurements, formed into a thin sheet using a cylindrical and focusing lens in 
combination for planar LIF measurements (PLIF), or used to illuminate a volume for three- 
dimensional or stereoscopic imaging. Reviews of the laser-induced fluorescence measurement 
technique are available from Eckbreth 42 and others . 43 ’ 44 The measurement technique works by 
inducing a transition, usually of an electron, from a lower energy state ( Ej ) to an excited energy 
state (Ei) via stimulated absorption of one or more photons in the atomic or molecular species of 
interest. In a two-level model assumption, the atom or molecule of interest in the fL state then 
returns to the Ej state by transferring energy via spontaneous emission of a photon 
(fluorescence), or by transferring energy non-radiatively through a collision with another atom or 
molecule (collisional quenching). 

For an atom, such as N or O, the energy required to induce an absorption transition of an 
electron to the instate from the ground (Ej) state is equal to the energy difference between the 
atom with an electron occupying the excited 
electronic orbital and the atom’s ground 
electronic orbital configuration, respectively. 

For molecules, such as N2 and O2, the energies 
of the E2 and Ej states include the energies 
associated with the vibrational and rotational 
motion of the molecule in addition to the 
energy associated with the molecule’s 
electronic configuration. Figure 2.1 shows a 
generalized two-level energy model for 



Figure 2.1: Two-level energy model of single- 
photon fluorescence. 



fluorescence with the stimulated absorption 
transition induced by a single photon. 

In Fig. 2.1, the rate at which the absorbing 
medium in the Ej state transitions to the Ei 
state is proportional to the product of the 
Einstein coefficient for stimulated absorption, 

B 12 , and the laser’s spectral irradiance (power 
per unit area per unit frequency), I v — IL V , 
where / is the irradiance (power per unit area) 
and L v is the laser spectral profile or laser line- 
shape (per unit of frequency). The energy of the absorbed and emitted photons is the frequency, 
v, multiplied by Planck’s constant, h. The Einstein A 21 and Bn coefficients describe the 
probabilities for emission and absorption, 45 respectively, while Q 21 is the collisional quenching 
rate. For absorbing species, the line-shape function, Y v , describes the spectral width for a 
particular energy level. This line-shape function combines broadening effects due to Gaussian- 
shaped Doppler broadening and Lorentzian-shaped homogeneous broadening mechanisms and is 
represented as: 42 


Figure 2.2: Two-level energy model of two- 
photon fluorescence. 


Y «=ikW na ' x} ai) 

Here, Av 0 is the Doppler-broadened line width due to the thermal motion of the absorbing 
species and V(a,x) is the Voigt integral function describing the convolution of homogeneous 
and Doppler broadening mechanisms. The line-width resulting from Doppler broadening is 

1 42 45 

given by. 
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Here, v 0 is the transition center frequency, c is the speed of light, k B is the Boltzmann constant, 
T the temperature, and m the species mass. The Voigt integral function is given by: 42 


V (a, x) = -J_° 


00 a 2 +(x-y ) 2 


dy 


(2.3) 


where a = (Av H /Av D ) A /ln(2) accounts for effects of homogeneous broadening ( Av H ) and 
x — {(v — v 0 )/ Av d }tJ 4 ln(2). Methods for computing this integral are given in Refs. 46-48. 


The integral of the product of the absorption line-shape function, Y v , and the laser’s spectral 
profile, L v , is defined by the overlap integral, G — f L v Y v dv, and describes what portion of a 
particular absorption transition is affected by the incident laser radiation. The rate constant for 


stimulated absorption, W 22 , which describes the rate at which species in Ej transition to E 2 via 
absorption of a single photon 42,44 is given by: 

W 12 = B 1Z IG (2.4) 

As the laser passes through a flowfield, it is absorbed at a rate corresponding to Eq. 2.4, 
inducing a transition between the Ej and E 2 states. Consequently, as the laser continues to 
propagate through the flowfield, the irradiance is continually diminished as a result of the 
absorption process. The Beer-Lambert law describes the magnitude of the decrease of the 
spectral irradiance, and is of the form: 49 


— dl v — 7 v /c s (v)dx 


(2.5) 


where k s is the absorption coefficient. The absorption coefficient can be related to B 12 by the 
relation / k s (y)dv — N S B 12 hv 0 /c , 49 where N s is the number density of the species, h is Plank’s 
constant, and c the speed of light. Integrating Eq. 2.5 results in the relation: 49 
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where I v 0 is the incident spectral irradiance and I v x the spectral irradiance after the beam travels 
a distance x. 


In certain circumstances, the energy separation between the Ej and E 2 states for the species 
being examined with LIF is such that the frequency of a photon necessary to induce an 
absorption transition is far into the ultra-violet portion of the electro-magnetic spectrum. 
However, achieving frequencies far into the UV with conventional laser systems can be 
difficult. 45 Additionally, strong absorption of UV radiation by many materials and gases 
(including air) below approximately 200 nm limits the application of laser systems capable of 
producing such frequencies. 42 One solution is to use two-photon LIF techniques to probe species 
having absorption transitions in the deep UV. With a two-photon technique, the frequency (and 
therefore energy) of each photon is half that required by a single photon for the same transition. 
Figure 2.2 shows a generalized two-photon fluorescence energy model for two-photon LIF. For 
atomic species, such as N or O, spontaneous emission is often observed between E 2 and an 
intermediate energy state, E 3 . 

For the two-photon LIF process in Fig. 2.2, the two-photon absorption cross-section results in 
a rate constant different from that in Eq. 2.4 and has the fonn: 42,45, 50-53 

W 12 = C 12 I 2 G (2.7) 

where Cn relates to the two-photon absorption cross-section. Note that the dependence of the 
rate constant scales with the square of the incident laser irradiance, I. This behavior arises from 
the probability of observing the simultaneous arrival of two photons, which is the square of the 
probability for the arrival of a single photon 45 


In both Figures 1 and 2, Q 21 is the so called quenching rate constant. It describes the rate at 
which energy is transferred through non-radiative collisions between excited atoms or molecules 
in the E 2 state and atomic or molecular collision partners of species i. This rate constant is 
computed similarly to that in Ref. 42 as: 

Q 21 = N T 'ZiXi°s,iVs,i (2.8) 

where N T — P/k B T is the total population of the excited state, x t is the mole fraction of 
quenching species i, o si is the collision cross-section between the excited species, 5 , and 

quenching species, i, and v s i — J 8k fJ T trans / ng s ( is the mean relative velocity 54,55 between the 

excited species and quenching species, with p s ( being their reduced mass and T trans their 
translational temperature. The spontaneous emission coefficient, A, in Figures 1 and 2 is also 
kn own as the Einstein A coefficient and describes the probability for spontaneous emission of a 
photon by an atom or molecule in the excited state. 45 

For the two-level model in Fig. 2.1, the rate of change of the populations N/ and N 2 , in the Ej 
and E 2 states, respectively, can be obtained from relations similar to those presented in Ref. 42: 

^ = — Ai fW 12 + N 2 (Q 21 +A 21 ) (2.9) 


JVi + JV 2 = JV 5 =* s / b JV t (2.10) 

Equation 2.10 is a conservation law saying that the combined populations of the Ej and E2 states 
is equal to the initial population of the excitation species’ Ej state, Ns, which is in turn the 
product of the species mole fraction, x s , the total population, N T , and the temperature-dependent 
Boltzmann fraction, f B . Equation 2.9 assumes that the laser intensity is sufficiently weak such 
that stimulated emission (W21) can be neglected. 


The Boltzmann fraction, f B , describes the ratio of the number of absorbers initially 
occupying Ej relative to all possible energy states at a particular temperature, T, when a system is 
in thermodynamic equilibrium. 55 This term has a general fonn given by: 42,55,56 



g } e-*i ,k * T 

I.jBj e~ E i ,kBT 


( 2 . 11 ) 


where g. is the degeneracy of state j. Degeneracy refers to the number of quantized states that 

exist in a given energy level, Ej . 55 The summation in the denominator is termed the partition 
function, Z, and for a diatomic molecule can be represented as the product of the individual 
partition functions for rotational and vibrational energies, Z = Z rot Z vib . The individual partition 
functions take a form similar to that in Ref. 55: 


( 2 . 12 ) 


Zrot — S;=0.5,1.5,2.5,- (2/ + l)exp [- k B T ro J 
Zvib = Zy=0,l,2,- exp 

where J is the rotational quantum number, v is the vibrational quantum number, Fj the rotational 
energy, G v the vibrational energy, T mt the characteristic rotational temperature, and T vi b the 
characteristic vibrational temperature. 

More involved differential rate equations can be formulated for the two-photon method like 
those presented in Ref. 53. It should be noted that the models presented in Figures 1 and 2 
neglect transitions from Ej to E 2 resulting from collisions, with a rate constant of Q12. This is 
usually a good assumption for large energy separations associated with the visible or UV 
transitions typically employed for LIF. Transitions from the E 2 and E 3 states resulting from pre- 
dissociation (Qpre) and ionization ( Q ion ), which are described in Refs. 42 and 44, have also been 
neglected. 


2.1 Linear, Steady State Solution 


If a continuous laser source is used to populate the E 2 state and detection of fluorescence 
occurs well after this source is turned on, then N 2 can be assumed to have reached its steady-state 
value. From this assumption, the left-hand-side of Eq. 2.9 is set to zero, resulting in two 
algebraic equations (Eqs. 2.9 and 2.10) for two unknowns (Nj and AT). The same assumption can 
be made to determine the population, AT, achieved by a pulsed laser source if the time required to 
reach steady state is short compared to the duration of the pulse. This population is computed as: 


N 2 


Xsf B N T W U 
W 12+Q 2 1 +A 21 


(2.13) 


The product of this excited state population and the spontaneous emission rate constant is 
N 2 A 2 i and represents the number of transitions per unit time and per unit volume. Integrating this 
constant value with respect to time gives the total number of transitions per unit volume during 
the detection period. Substituting Eq. 2.4 for W 12 , and assuming that Wj 2 is small compared to 
A 21 and Q 2 i (which is valid for low-intensity excitation), a relation for the total number of 
photons collected via fluorescence by the detection device is : 44 

Slip = Xsf B N T B 12 IG ( pt det V (2.14) 

where <t> = A 21 /(Q 21 + A 21 ) is the fluorescence yield, t de t is the period of detection, V is the 
volume probed by the laser source, fl is the solid angle over which detection occurs, and 77 is the 
detection efficiency. The fluorescence yield, O, describes the fraction of de-excitation transitions 
that occur via spontaneous emission (fluorescence) relative to all de-excitation transitions (i.e. 
spontaneous emission and collisional quenching, assuming pre-dissociation and ionization are 
negligible). Neglecting the constants in Eq. 2.14, a generalized form for the fluorescence signal 


similar to that presented in Refs. 44 and 57, including its thermodynamic (x,P, T) dependencies 
and velocity dependence (arising from the Doppler effect), is: 

Slip « Xsf B (T)B 12 IG( Xs > P. T, U)<$>(x s ,Xi. P. T)t det (2.15) 


2.2 Non-Steady State Solution 

For a pulsed laser source, if N 2 is changing during the period of detection, then the entirety of 
Eq. 2.9 must be solved to obtain the time-dependent population, N 2 (t). Similarly, if the duration 
of the laser pulse is of the same order of magnitude as the time required to reach steady state, 
then the population at the end of the laser pulse, N 2 (t laser), can be obtained by solving Eq. 2.9. By 
assuming the laser intensity behaves as a Heaviside step function in time, the population N 2 (t), 
existing during laser excitation can be calculated by substituting the relation for Ni from Eq. 2.10 
into Eq. 2.9 and integrating with respect to time: 44 

N 2 (t) = Xsf B N T ^-( 1 - e~ rt ), 0 <t< t laser (2.16) 

where r = W 12 + Q 2 \ + d 21 . The inverse of this value, r, is the characteristic time needed to 
achieve steady state. This solution assumes an initial condition which typically specifies the 
initial excited state population to be zero (A?(0) = 0). 

When the laser source is turned off, the only pathways to return to the ground state from the 
excited state are assumed to be through either spontaneous emission or collisional quenching. 
Therefore, for the period following laser excitation, the first term on the right-hand-side of Eq. 
2.9 is zero. This modified version of Eq. 2.9 is then used to determine the excited-state 
population for the period after laser excitation by integrating with respect to time: 44 

N 2 (t) = Xsf B N T ■ (1 - eraser) ■ e -P-Gaser) frllF f t > t laser (2.17) 


This solution assumes that the laser intensity is turned off instantaneously, with an initial 
condition given by Eq. 2.16, evaluated at t = t iaser . This solution shows that the population in 
the excited state decreases exponentially in time after the laser pulse. 

The term t lif — (Q 2 \ + d 2 i) _1 i n Eq. 2.17 is referred to as the fluorescence lifetime and 
describes the rate at which the population in a particular excited state transitions to a lower state. 
As with the derivation of Eq. 2.14, the total number of photons collected via fluorescence over 
the detection period is determined by multiplying Eqs. 2.16 and 2.17 by the spontaneous 
emission rate Einstein coefficient A 21 , accounting for the collection volume and detection 
system, and integrating with respect to time: 44 
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Like Eq. 2.14, this solution assumes that Wn is small compared to A 2 i and Q 21 . However, if 
the laser irradiance, /, is sufficiently large such that W 12 and W 2 j are of the same order of 
magnitude as Q 21 and A 2 i, then both must be included in the solution. This gives r — W 12 + 
W 21 + Q 21 + A 21 during laser excitation. By defining the irradiance at which the fluorescence 
signal saturates as 42 ’ 44 I sat — (Q 21 + A 21 )/(_B 12 + 5 2 i) an d rearranging r such that r — 
t lif\I /I sat + 1], a more detailed formulation for the fluorescence signal, including effects from 
laser saturation, is given as: 
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When / « I sat , Eq. 2.19 simplifies to a fonn similar to that of Eq. 2.16. 


2.3 Multi-Level Fluorescence Modeling 

In the preceding discussion, a simplified two-level model of fluorescence was put forth, 
which provides for an understanding of the most basic physical mechanisms and energy transfer 
processes involved. This basic description allows for the development of a general analytic 
relation between the fluorescence signal and these mechanisms and processes. However, such a 
simple description does not account for rotational or vibrational energy transfer occurring 
between the absorbing species and the surrounding gas mixture. Since each electronic energy 
level depicted in Fig. 2.1 can have vibrational and rotational fine structure for molecular species, 
such energy transfer processes result in a redistribution of the populations to various vibrational 
and rotational states. To account for these processes, analytical models including multiple energy 
states and a system of population rate equations are used. Discussion of such modeling aspects 
can be found in Refs. 42-44. 


2.4 Translational Temperature 

Translational temperature resulting from the random thermal motion of the absorbing species 
can be determined by measuring the Doppler-broadened width of the absorption line-shape 
function of Eq. 2.1. If the temperature, pressure, and mole fractions at a particular location are 
assumed to be constant, then the Boltzmann fraction and fluorescence yield in Eq. 2.15 at that 
location are also assumed to be constant. If these assumptions are valid, and if the laser line- 
shape function is known, then the fluorescence signal will be proportional to the overlap integral 
at a particular frequency. By scanning the laser’s spectral frequency across a particular transition, 
the absorption line-shape function can be obtained. 

Once the absorption line-shape is obtained, the contribution of homogeneous broadening 
mechanisms to the line-shape must be determined to allow for calculation of the Doppler- 
broadened line width. Homogeneous broadening can include components from Van der Waals 

CO 

broadening and Stark broadening effects, in addition to collisional pressure -broadening effects. 
In many cases, the Van der Waals and Stark effects can be ignored, resulting in homogeneous 


broadening from pressure alone. If the pressure is known, then the component of the line-shape 
due to homogeneous broadening can be determined. Alternately, if the pressure is sufficiently 
low, then collisional pressure broadening can be ignored. A fitting algorithm can then be used to 
determine the contribution of Doppler broadening to the line-shape. This allows for a more 
accurate temperature calculation as only one parameter, Av D , must be iterated upon to lit the 
line-shape. The translational temperature, T tram , of the absorbing species can then be obtained by 
modifying Eq. 2.2: 


1 trans 


m (cAv d > 
8ln(2 )k B V v 0 > 


( 2 . 20 ) 


Both single-photon 58 ' 60 and two- 
photon 61 ' 6 '^ methods have been used to 
measure translational temperature. The 
requirement of scanning a laser’s 
frequency to resolve the Doppler- 
broadened line-shape can limit the time 
resolution of the temperature 
measurement such that only average 
translation temperatures can be obtained. 


Figure 2.3 shows LIF data (dotted 
line, upper plot) from Ref. 58 obtained 
by performing a frequency scan with a 
semiconductor laser source over an argon 
metastable transition centered at 810.4 
nm. Each frequency scan over the 
absorption transition took 500 seconds to 
complete. In the experiment, argon 
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Figure 2.3: Experimental LIF data measured (dotted 
line, upper plot) in an Argon plasma and corresponding 
fit of Voigt line -shape (solid line, upper plot) to the 
experimental data (upper plot) and difference between 
experimental data and fit (lower plot). Image reprinted 
from Ref. 58 with permission of the authors and the 
publisher. 

plasma was generated using an inductively-coupled plasma torch. The translation temperature, T , 
and electron number density, n e , were obtained by fitting the Voigt line-shape function, which 
accounted for both Doppler-broadening and several collisional broadening mechanisms, to the 
experimental LIF data (solid line, upper plot). The lower plot in Fig. 2.3 shows the residual: the 
difference between the experimental LIF data and the fitted line-shape. An analysis of one of the 
scans in Ref. 58 gave T tran s = 6800±690 K, resulting in a -10% measurement uncertainty. 


2.5 Rotational Temperature 

For molecules, the rotational temperature, T rot , can be obtained by exciting two absorption 
transitions and relating the signal intensities observed from each transition to the rotational 
temperature through modification of Eq. 2.15. This method was demonstrated in Ref. 66 using a 
two-photon excitation method to measure rotational temperature at a point in a cold (-300 K) 


turbulent flow. For each rotational transition, the population fraction, f B , is dependent on the 
rotational energy, Fj, rotational quantum number, J, and rotational temperature, T rot , of the 
absorbing species. This assumes that excitation occurs in the same vibrational state, v, such that 
the vibrational energy, G v , is constant for each probed rotational transition. If/ S , G, O, and t det 
are assumed to be independent of the rotational state of the absorbing species, then the ratio, R , 
of fluorescence signals can be simplified to the following relation, similar to that in Refs. 44, 67- 
70: 


s LIF,i _ r _ q b 12 ,i E i (2/i+l)exp[ Fj i/ksTrot] 
s LIF,j B i 2,jE j (2/; + l)exp[ F j j /kBT ro t\ 


( 2 . 21 ) 


where the subscripts i and j refer to the particular excited absorption transition, E is the laser 
energy such that E — ati aser I, a is the beam cross-sectional area, ti aser is the temporal pulse 
width of the laser, and C is a constant which includes all terms that are independent of the 
rotational state of the absorbing species. Solving for T ro , in Eq. 2.21 yields: 
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Here, A E rot — ( Fj :i — Fjj ), and is the difference in energy between the probed rotational levels. 
To minimize the uncertainty associated with the measured rotational temperature, 8T rot , a 
propagation-of-error analysis can be performed on Eq. 2.21 by computing the derivative of the 

z:o n i 

signal ratio, R, with respect to T rot : ’ 
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This analysis suggests that by increasing A E rot , and thereby probing two widely separated 
rotational energy levels, the uncertainty 
in T mt can be proportionally decreased 
for a given error in R. 68 

Figure 2.4 shows a mean rotational 

temperature map of nitric oxide 

measured in the wake of a 6. 3 5 -nun- 

thick flat plate in a Mach 3 turbulent 

flow, taken from Ref. 71. The map was 

obtained by averaging a series of images 

taken for excitation from rotational 

levels J = 8.5 and J = 10.5, computing , ■ , ■ , 

figure 2.4: Rotational temperature map obtained m the 

the ratio of the two averaged images, and supersonic wake of a flat plate using two-line rotational 

using a relation similar to that in Eq. thermometry. Image reprinted from Ref. 71 with 

2.22 to compute T mt . In this experiment, permission from the authors and Springer Science and 

Business Media. 




it was shown that the turbulent nature of the flowficld required excitation of two relatively 
closely spaced rotational levels (J = 8.5 and 10.5). This requirement arose from the nonlinear 
relation between R and T rot in Eq. 2.22, which can be heavily influenced by flowfield turbulence. 

Figure 2.4 shows a time-averaged temperature measurement obtained using a single laser, 
which probed the two transitions in separate wind tunnel runs. To obtain instantaneous 
temperature measurements in a plane using this method, a two-laser, two-camera system can be 
used, as in Refs. 67-69. In Ref. 68, instantaneous two-line temperature imaging of OH in a shock 
tube at a nominal pressure of 40.53 kPa and temperature range of 1500-2950 K resulted in 
temperature measurements with errors of ~20%-25%. 

Multi-rotational-line temperature methods also exist and can potentially provide relatively 
higher sensitivity and dynamic range T rot measurements, particularly at lower temperatures 
where rotational energy levels are closely spaced. " Such methods are appropriate when large 
variations in T rot are expected, requiring probing of several rotational levels. The rotational 
temperature is computed by exciting multiple rotational lines and measuring the signal. This 
measured fluorescence signal, together with the measured laser energy are then substituted into 
the following relation: 70 
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A plot of this logarithmic term versus Fjj, kn own as a Boltzmann plot, yields a linear trend with 
a slope of — k B Trot- Refs. 72-75 used this method to characterize rotational temperatures in arc 
jet flows. An alternative multi-line approach used in Ref. 76 for measurements in a flame 
involved fitting a simulated rotational spectrum to experimental data to obtain T rot . 


2.6 Vibrational Temperature 

The vibrational temperature, T vi b, can be measured using an approach similar to that used for 
rotational temperature. For T vi b measurements, a two-line approach can be used. Such an 
approach is presented in Ref. 77, in which the same rotational levels, J t — Jj, are probed in two 
different vibrational levels, A Vj , resulting in nearly constant rotational energies, Fj, but 
different vibrational energies, G v . Assuming Xs-> G > , and t det to be independent of the 

rotational and vibrational levels, the ratio of fluorescence signals, R, can be used to measure 
vibrational temperature, T vi b, in a way analogous to that in Eqs. 2.21 and 2.22. 

As with rotational temperature measurements, a multi-vibrational-line temperature 
measurement method can also be used to infer T vi b by modifying the logarithmic term from Eq. 
2. 24: 70 
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If the rotational level, / p is fixed and the rotational energy, Fj t , is kept approximately constant, 
then the vibrational temperature can be computed independently of rotational temperature. This 
is done by generating a Boltzmann plot of vibrational energy, G v i , versus this logarithmic term, 
resulting in a linear relation with a slope equivalent to — l/{k B T vib ). 


Figure 2.5(a), taken from Ref. 70, shows a Boltzmann plot generated using a relation like that 
in Eq. 2.25 to infer vibrational temperature. The image data in the experiment was used to 
measure T vi b at a point downstream of a free-piston shock tunnel nozzle on the centerline of the 
flow. In this experiment, three absorption transitions were selected that had nearly constant 
rotational energies, Fj, but differing vibrational energies, G v . In Fig. 5(b), also taken from Ref. 
70, the vibrational temperature of NO was observed to be constant with distance downstream of 
the nozzle, even though the nozzle was conical and the rotational temperature was shown to be 
decreasing with distance downstream, as the flow expanded. A simple computational model 
(shown in Fig. 5(b)) was not able to accurately predict the measured vibrational temperature. 70 
As shown in Fig. 5(b), T vlh = 785±30 K, giving an uncertainty of ~4%. 
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Figure 2.5: Vibrational temperature measurement obtained by probing multiple vibrational levels and 
generating (a) a Boltzmann plot and (b) vibrational temperatures plotted as a function of distance 
downstream. Images taken from Ref. 70 with permission of the authors. 


2.7 Species Concentration 

Returning to the simplified Eq. 2.14, it is possible to quantify the number density, ^ s iV r , 
using LIF. However, several factors complicate the interpretation of a LIF signal and its relation 
to concentration. For instance, Q 21 in Eq. 2.8 depends on the concentration and temperature of 
the absorbing species and all other constituent species in the probed mixture. If the temperature 
of the mixture were unknown, then a separate measurement of temperature would be required to 
determine both Q 21 and the temperature-dependent Boltzmann fraction. The measurement system 
also requires an appropriate absolute intensity calibration to a known reference, which can be 
further complicated by the collection optics and filters used in such a measurement. The response 


of the detection system as a function of S LIF intensity and spectral frequency must also be well 
understood. The following discussion highlights some of the methods that have been used to 
provide concentration measurements that address some of these issues. 

2.7.1 Saturated LIF 

To avoid the complications associated with an unknown quenching rate constant, a relatively 
high laser irradiance where / » I sat can be used to probe the absorbing species . 42 In this 
instance, the stimulated absorption and emission rates dominate both collisional quenching and 
spontaneous emission rates; that is, (W 12 + W 21 ) » ( Q 21 + A 21 ). The result of this leads to a 
modified fonn of Eq. 2.14 where stimulated emission is included and the fluorescence signal is 
of a fonn similar to that in Ref. 42: 

Slip = Xsf B N T A 21 -^-t det V £17 (2.26) 

# 12+^21 471 

If A 21 , B 12 , and B 21 are known, then this method can be used to measure concentration if the 
remaining terms, such as f B , can be determined. In practice, it is difficult to achieve / » I sat 
because high powered pulsed laser beams are typically Gaussian, both spatially and temporally, 
so much of the acquired signal can originate from the lower intensity edges of the Gaussian 
beam, which excite the fluorescence linearly. Further discussion of this method is found in 
Refs. 42 and 43. 

2.7.2 Pre-dissociation LIF 

Pre-dissociative fluorescence involves inducing a transition to a level in the excited state, 
which can then couple to a dissociative state, resulting in dissociation of the molecule into 
smaller molecules or atoms. Pre-dissociative fluorescence concentration measurements are 
similar to saturated fluorescence measurements in that the dependence on Q 21 in Eq. 2.14 can be 
removed, in this case by exciting the absorbing species to a pre-dissociative state where the pre- 
dissociation rate, Q pre , is faster than both Q 21 and A 21 . A discussion of this mechanism can be 
found in Refs. 78 and 79 and a general discussion of the method can be found in Refs. 42 and 43. 
An application of the pre-dissociative technique was presented in Ref. 80 to measure time- 
averaged OH concentrations in a supersonic hydrogen-air turbulent combusting flow. 
Concentrations as high as 5.4x10 cm" with uncertainties of 21% or less were reported. 

2.7.3 Short-Pulse LIF 

There are two methods of using short laser pulses to quantitatively determine concentration. 
The first probes the species of interest with a laser pulse having a duration much shorter than the 
time between collisions with other species. With this method, the concentration can be 
determined by relating the measured exponential fluorescence decay behavior to the initial signal 
magnitude during the short-pulse excitation period. While this decay is a function of both A 21 and 
Q 21 , extrapolation of the exponential behavior to initial excitation provides an inferred peak 


intensity magnitude, which is assumed to be independent of collisional quenching effects. 4. This 
independence is assumed since a sufficient number of collisions, required to transfer energy non- 
radiatively from the excited state, would not yet have occurred. A discussion of this method can 
be found in Ref. 43. Reference 81 uses a form of this short-pulse method for OH concentration 
measurements in a turbulent flame. 

A second method for measuring concentration involves using a laser pulse that is shorter than 
or comparable to the fluorescence lifetime but not shorter than the collisional timescales 
described in Ref. 42. In this approach, the fluorescence lifetime is measured directly, usually 
using a photomultiplier tube. This short-pulse LIF method can also be used to determine 
spontaneous emission " and collisional quenching ’ rate constants when the thermodynamic 
conditions of the probed mixture are known. Even for u nkn own conditions, measurement of the 
fluorescence lifetime allows for the quenching rate to be determined and accounted for. Refs. 85 
and 86 use a two-photon LIF technique in which fluorescence lifetime measurements are used to 
correct for collisional quenching effects in a similar manner. Figure 2.6(a), taken from Ref. 85, 
shows spectral scans of atomic nitrogen in the NASA Ames Aerodynamic Heating Facility arc 
jet flow (red data) and laboratory flow reactor (green data). Fluorescence measurements from the 
flow reactor and a krypton reference cell were used for intensity calibration in order to provide 
absolute atomic nitrogen number densities. The fluorescence lifetimes observed in the arc jet and 
flow reactor were used to correct for collisional quenching effects. Figure 2.6(b), also taken from 
Ref. 85, shows radial atomic nitrogen number density measurements for two air/argon arc jet 
runs taken 35.6 cm downstream of the nozzle exit. The reported uncertainty in number density 
was ~12%. 85 




Figure 2.6: (a) Spectral scans of atomic nitrogen in the NASA Ames Aerodynamic Heating Facility arc 
jet flow and in a laboratory flow reactor and (b) measured number densities spanning the radial direction 
of the arc jet 35.6 cm downstream of the nozzle exit. Image taken from Ref. 85 with permission of the 
authors. 


A similar approach is being developed for mole fraction measurements in hypersonic 

0-7 OO 

turbulent boundary layers using naphthalene PLIF. ’ One benefit of using naphthalene is that it 
sublimates at slightly elevated temperatures (with respect to room temperature), allowing for the 
study of scalar transport effects in transitional boundary layers, turbulent boundary layers, and 
ablating surfaces. 

2.7.4 Bi-Directional Beam LIF 


If two overlapping, counter-propagating beams are tuned to the same transition, the ratio of 
fluorescence signals results in the cancellation of all terms in Eq. 2.15 at a point in the flow, with 
the exception of the spatially dependent laser irradiance. Thus, the spatially varying ratio of 
fluorescence signals can be equated to the spatially varying ratio of irradiances of the respective 
beams. The Beer-Lambert law of Eq. 2.5 can then be used to relate irradiance, and therefore 
fluorescence signals, to number density via the relation: 49 

k s (v) = N s (p s (y) (2.27) 

where (p s (v) is the transition cross-section of the absorbing species, describing the absorbing 
species’ interaction with the incident irradiance on a per atom or molecule basis. 45 Combining 
Eq. 2.27 with Eq. 2.6 allows for a relation between the natural logarithm of the signal ratio and 

on 

the absorbing species number density to be made: 
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where Slif, i (x) and Slip, 2 (x) are the spatially varying fluorescence signal intensities of the 
forward- and backward-propagating beams, respectively. In Eq. 2.28, cp s 0 is the peak transition 
cross-section, with cp s 0 — cp s (y 0 ). 45 


This concentration measurement technique was first demonstrated in Ref. 89. The technique 
only requires that a value for cp s 0 be known in order to make an absolute concentration 
measurement of N s . For total absolute concentration, however, f B must either be assumed or 
measured, as from Eq. 2.10, N s — Xs/b^t- Figure 2.7, taken from Ref. 89, graphically depicts 
the process of obtaining a concentration of OH in a methane-air flame. In the top plot, the 
fluorescence intensities are plotted versus position across the flame. The middle plot shows the 
natural logarithm of the ratio of these two signals with position and the bottom plot shows the 
calculated absolute OH concentration using Eq. 2.28. Measurements of OH concentration were 
also made in a plane in Ref. 89 by imaging fluorescence from two counter-propagating laser 
sheets. Equation 2.28 was applied using each row of pixels to compute signal ratios. Single-shot 
imaging of a hydro gen/air/oxygen flame allowed for calculation of means and standard 
deviations of OH concentration. Using 10 images for each laser sheet direction (20 total images), 
Figure 2.7 of Ref. 89 showed a mean peak OH concentration of ~8.8xl0 14 cm' 1 with a standard 
deviation of ~±1.7xl0 14 cm' 1 , which is — 19.7% of the mean. 


2.7.5 Combined LIF/Rayleigh/Raman 



In(Ratio) 



In reacting flows, flows undergoing 
transition-to-turbulence, and turbulent 
flows, the time -varying nature of flow 
properties makes it difficult to ascertain 
concentration. By making simultaneous 
time-resolved measurements of 
temperature and major species 
concentrations, the thermodynamic - 
dependent parameters in Eq. 2.15 
(Boltzmann fraction, overlap integral, and 
fluorescence yield) can be determined, and 
a direct relation between concentration 
and fluorescence signal established. A 
combined fluorescence, Rayleigh 
scattering, and Raman scattering 
measurement technique was used for this 
purpose in Ref. 90. In that experiment, the 
Raman scattering was used to measure 
instantaneous concentrations for major 
species (O 2 , N 2 , H 2 0, and H 2 ) in a 
turbulent hydrogen jet flame. This allowed 
for the mole fractions of the collision 
partners, 2 % in Eq. 2.8 to be calculated. 

The combined Rayleigh/Raman scattering 
measurements were then used to measure 
temperature. The temperature-dependent 
collisional cross-sections and relative 
velocity in Eq. 2.8 could then be 
calculated, and the collisional quenching 
rate, Q 2 i, determined. The temperature- 
dependent absorption line-shape and 
Boltzmann fraction could also be 
determined from the temperature 
measurement. The additional measurement of laser energy then provided for a direct relation 
between NO fluorescence signal and concentration to be established. In Ref. 91, this technique 
was used to make concentration measurements of CO in turbulent premixed and stratified 
CfU/air flames. A stated CO concentration accuracy and precision (single standard deviation) of 
10% and 4.5%, respectively, were given. 91 
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Figure 2.7: Graphical depiction of methodology used 
to measure concentration using bi-directional beam 
LIF. Top plot shows spatially varying fluorescence 
signals, middle plot shows natural logarithm of ratio of 
these signals, and bottom plot shows OH concentration 
calculated using Eq. 2.28. Image reprinted from Ref. 
89 with permission of the authors and Springer Science 
and Business Media. 


2.8 Doppler-Based Velocimetry 


The Doppler Effect can be used to determine flow velocities for atomic and molecular 
species. The translational motion of the absorbing species in the direction of the excitation 
laser’s propagation, described by a velocity component U, results in a shift of the absorption 
line-shape function away from its transition center frequency, v 0 , according to: 49 

Av = % 0 (2.29) 

This velocity-dependent frequency shift of the absorption line-shape function is implicit in 
the overlap integral, G, and thus, its effect on fluorescence signal can be seen through the 
dependence of G on the velocity component, U, in Eq. 2.15. When a component of translational 
motion of the absorbing species opposes the laser’s direction of propagation, corresponding to 
the negative (— ) solution of Eq. 2.29, the incident laser radiation appears to be at a higher 
frequency from the perspective of the gas. Hence, as the laser’s frequency is scanned over the 
absorption transition, Y v , the measured intensity of fluorescence corresponding to this profile is 
shifted toward a lower frequency, or red-shifted. The converse is true when motion is in the same 
direction as the laser’s propagation, corresponding to the positive (+) solution of Eq. 2.29, where 
the absorption profile is shifted toward a higher frequency, or blue-shifted. 

In one implementation, velocity component measurements can be made by scanning the laser 
over an absorption transition in both the measurement volume and a reference cell. In this case, 
the Doppler shift between the absorption profiles is used to compute an average velocity 
according to Eq. 2.29. Examples of such measurements include those in an arcjet, 59,60,64,92 
supersonic underexpanded jets, 93 ' 97 shock tunnel, 98,99 and non-reacting supersonic flow with a 
rearward-facing step. 100 ’ 102 Such measurements require that the flowfield be relatively steady 

since shot-to-shot fluctuations in 
fluorescence intensity from 

thermodynamic and fluid dynamic 
variations can affect the measured line- 
shape. The velocity measurements can 
also be complicated by a frequency shift 
in the line-shape function resulting from 
collisional effects (pressure). In 
compressible flows, the pressure field 
can vary significantly, and hence the 
collisional shift in the line-shape function 
can likewise vary. Absorption of laser 
energy can also result in an apparent 
frequency shift in the line-shape 
function. Recall that the line-shape 
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Figure 2.8: Two-component velocity measurement 
using fluorescence-based Doppler velocimetry 
technique. Vectors denote flow direction and color map 
represents magnitude. Image reprinted from Ref. 99 with 
permission of the authors and Springer Science and 
Business Media. 


function must be inferred from the overlap integral, G, which represents the convolution of the 
absorption line-shape and laser spectral line-shape. If absorption is significant, then irradiance 
will vary spatially according to Eq. 2.6, resulting in a spatially varying overlap integral. An 
analysis of how absorption affects the line-shape frequency shift is presented in Ref. 99. If flow 
symmetry is assumed, then the frequency shift due to pressure can be estimated. Flow 
symmetry can also be used in the application of Eq. 2.6 to correct for the frequency shift 
resulting from absorption. 60 Alternatively, if two counter-propagating laser beams are used, these 
frequency shift effects can be completely removed. The use of counter-propagating beams results 
in two excitation peaks, separated in frequency by twice the velocity-induced Doppler shift. The 
need for a reference cell measurement to ascertain velocity can also be removed, as a counter- 
propagating beam approach is self-referencing. 95 ’ 103,104 

Figure 2.8 in Ref. 99 shows a two-component velocity measurement about a heat shield 
model taken in a hypersonic free-piston shock tunnel. In this experiment, laser sheets were 
directed in both the radial (vertical) and axial (horizontal) directions so that measurements of the 
Doppler-shifted absorption profdes for the respective directions could be obtained and compared 
with measurements from a static reference cell to infer velocities. An estimation of velocity 
errors incurred from the frequency shifts due to absorption and collisional effects was also 
perfonned. In Ref. 99, three separate absorption transitions were probed to measure two different 
velocity components, with one transition being used for both components. This resulted in a 
measured axial freestream velocity of 2394±68 m/s and a measured radial velocity of 53±50 m/s, 
giving respective uncertainties of -2.8% and -94.3%. 99 

Another form of the fluorescence-based Doppler velocity measurement is a fixed frequency 
method, which can allow for an instantaneous velocity component measurement. With the fixed 
frequency technique, a narrow linewidth laser is tuned off the absorption profile peak to a point 
where the slope of the profile is maximum, as described in Refs. 105 and 106. Assuming that the 
absorption profile is approximately linear in the region of maximum slope, the measured signal 
intensity can be related to the Doppler shift of the profile. This fixed frequency Doppler 
velocimetry technique has been applied to a free jet, supersonic underexpanded jet, and 
reacting supersonic flow. 109 The stated random and systematic errors in Ref. 108, when added in 
quadrature, gave a total uncertainty of -12%. In Ref. 109, the stated lowest time-averaged and 
single-shot uncertainties achieved were -3% and -15%, respectively, for a 1600 m/s velocity 
range. 

2.9 Flow-Tagging Velocimetry 

Another technique by which velocity can be measured using fluorescence is flow-tagging 
velocimetry. Fluorescence-based flow-tagging velocimetry is a time-of-flight technique that 
involves laser excitation — or tagging — of the gas along a line, series of lines, or grid pattern. 
With this form of velocimetry, the species of interest in the gas absorbs the incident radiation 
from a laser source, which induces either of the following: 1) fluorescence, 2) a reaction that 


forms a product that then emits a photon via fluorescence, or 3) a reaction that forms a product 
that can then be probed with another laser source to induce fluorescence. Images of the 
fluorescence pattern are acquired at two time delays, with velocity computed by measuring the 
displacement of the tagged molecules between images. Typically, a line or series of lines can be 
used to measure a single-component of velocity while a crossed grid pattern can be used to 
measure two-components. Two general fluorescence-based methods of flow-tagging velocimetry 
are discussed here; one that requires a single laser source and another that requires multiple laser 
sources. The main advantage of flow-tagging velocimetry, as compared to most Doppler-based 
methods (which are time averaged), is that it can make instantaneous (single-shot) measurements 
with fast time resolutions (as short as a few hundred nanoseconds). A disadvantage, however, is 
that flow-tagging velocimetry cannot provide full velocity field information. A broader 
discussion of molecular-tagging velocimetry, which relies on molecular tracers for flow-tagging, 
is provided in Ref. 110. 

2.9.1 Single-laser methods 

The first method involves either direct or indirect excitation of fluorescence with a single 
laser source. The first application of this method to a gaseous flow involved excitation of 
phosphorescence of biacetyl molecules, as described in Ref. 111. If this method is used, the 
fluorescence lifetime of the tagged molecules must be long enough so that advection provides for 
measurable displacements with tagged regions having signal intensities that are still above the 
detection limit of the imaging system at the time the second image is acquired. Typical 
experiments involve capturing a reference image acquired during the tagging process, or a 
relatively short time thereafter. If a single-framing camera is used, such as in Ref. 112, a single 
reference image or set of reference images is acquired. The timing of the single-frame camera is 
then delayed and a subsequent image or 
set of images is then acquired. If a two- 
camera system or dual-framing camera is 
used, such as in Ref. 113, the delayed 
image is acquired in sequence after the 
reference image. The velocity is computed 
by measuring the displacement of the 
tagged molecules that occurs in the time 
between when the reference and delayed 
images were acquired. The form of this 
technique relying on direct excitation of 
fluorescence for flow tagging has been 
applied to the study of supersonic jets, 114 ' 

120 hypersonic boundary layers, 112 ’ 113 ’ 121 ’ 122 
and arcjet flowfields. 60 



Figure 2.9: Single-line excitation of nitric oxide 
fluorescence used to study hypersonic boundary layer 
flow over a flat plate. Images, from left to right, 
correspond to camera delay settings of 0 ns, 250 ns, 
500 ns, and 750 ns. Image reprinted from Ref. 112 
with permission of the authors. 


Figure 2.9, taken from Ref. 112, shows images from tagging a single line of nitric oxide 
using direct excitation of fluorescence within a flat plate hypersonic laminar boundary layer. In 
this Figure, the left-most image corresponds to the reference image, while the remaining images, 
from left to right, correspond to delayed images taken at 250 ns, 500 ns, and 750 ns after tagging, 
respectively. Measurements of freestream velocity, spatially averaged from a point just above the 
velocity boundary layer (3 mm) to 15 mm above the flat plate, resulted in a mean of 3,035±100 
m/s at 90% confidence, giving an uncertainty of 3.3% of the mean. Single-shot uncertainty 
estimates for a 3,000 m/s freestream flow and for camera delay settings of 250 ns, 500 ns, and 
750 ns were 4.6%, 3.5%, and 3.5%, respectively. 112 

An indirect excitation scheme, as described in Refs. 123-125, relies on photodissociation of 
molecular nitrogen for flow tagging. The technique uses a femtosecond laser pulse to dissociate 
molecular nitrogen into two nitrogen atoms, which then recombine after a collision, forming 
molecular nitrogen in an intermediate state. A subsequent collision brings the molecular nitrogen 
to an excited electronic B state, which then emits a photon via fluorescence upon transitioning to 
the excited electronic A state. Ref. 125 provides a description of this process. One benefit of this 
indirect technique, kn own as Femtosecond Laser Electronic Excitation Tagging (FLEET), is that 
the recombination rate of dissociated atomic nitrogen allows for a much longer fluorescence 
lifetime. This would allow displacements to be measured over greater time scales, providing for 
accurate measurements of velocity in low-speed flow regions, such as in a hypersonic wake flow. 
Additionally, the technique relies on molecular nitrogen for tagging, which is present in most 
hypersonic facilities. 

2.9.2 Multi-laser methods 

A second method of fluorescence-based flow-tagging velocimetry involves writing a line, 
series of lines, or grid pattern into the flowfield by one of several laser-based mechanisms. This 
pattern can then be interrogated, or read, by subsequent laser pulses to induce fluorescence, 
allowing for the determination of velocity through measurement of the displacement of the 
pattern. Such techniques usually involve two or three different lasers and are therefore more time 
consuming to set up and more difficult to execute. 

One mechanism by which a pattern can be written into the flowfield is via ionization of the 
absorbing species, kn own as Laser Enhanced Ionization (LEI) flow tagging. The tagging process 
is accomplished by promoting the species (such as sodium in Refs. 126-128) to a higher energy 
state via laser excitation near the ionization limit. Collisions then result in the ionization of the 
species, with the tagging pattern corresponding to the ionized regions. Subsequent laser pulses 
are used to induce fluorescence of the absorbing species in regions that have not been photo- 
ionized. Supersonic measurements of velocity in a shock tube were performed using LEI flow- 
tagging in Refs. 126 and 127, and hypersonic velocity measurements in an expansion tube were 
perfonned in Ref. 128. 


A second mechanism that can be used to write a pattern into the flowficld is by vibrational 
excitation of molecular oxygen via Raman pumping. This is followed by reading the pattern of 
the vibrationally-excited oxygen by inducing fluorescence. The technique, known as Raman 
Excitation and Laser-Induced Electronic Fluorescence (RELIEF), is advantageous as it relies 
on the flow-tagging of oxygen which is a common working gas of most hypersonic facilities. 
This technique has been used to characterize turbulence in a free jet and underexpanded jet. 
The RELIEF technique, however, is limited to temperatures below 750 K. Above this 
temperature, a significant fraction of oxygen molecules are vibrationally excited, making it 
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difficult to distinguish the tagged molecules from the background. 

Yet another mechanism involves using one laser to photo-dissociate a molecular species. 
This results in the formation of a product species for which a second laser can be used to read the 
location of the written pattern by exciting laser-induced fluorescence in the product species. The 
formation of the product species typically occurs through one or more reactions. A list of partner 
species used in the writing and reading process include: H 2 0-0H, 133-139 N 2 0-N0, 14 ° 0 2 -03, 141,142 
N 2 /0 2 -N0, 143-146 and N0 2 -N0. 116,147-154 

Figures 10(a) and 10(b), taken from Ref. 151, show images of vibrationally excited NO 
fluorescence in a grid pattern fonned via photodissociation of N0 2 in a mixture containing 6.3% 
N0 2 in N 2 . The images were taken in a supersonic underexpanded jet 400 ns (Fig. 2.10(a)) and 
800 ns (Fig. 2.10b) after the pattern was written into the flow with a two-dimensional array of 
355 mn beams. Two components of velocity were obtained by relating the displacement of the 
grid in the left image to a grid imaged in a stationary gas. The upper half of Fig. 2.10(c), taken 
from Ref. 151, shows measured streamwise velocities compared with computation, shown in the 
lower half of Fig. 2.10(c). The use of two pulsed dye lasers in this experiment pennitted the 
excitation of both a low (/ = 1.5) and high (/ = 8.5) rotational level within the same vibrationally 
excited (v = 1) NO state. The fluorescence images corresponding to the low-/ (Fig. 2.10(a)) and 
high-/ (Fig. 2.10(b)) rotational levels also allowed for measurement of the rotational temperature 
by using a calculation similar to that presented in Eq. 2.22. The top half of Fig. 2.10(d), from 


(a) (b) (c) (d) 

Figure 2.10: Flow-tagging images of vibrationally excited NO fluorescence obtained (a) 400 ns and (b) 
800 ns after photodissociation of N0 2 . The two-dimensional grid pattern allows for calculation of two 
velocity components. Measured (c, top) and computed (c, bottom) streamwise velocity maps. Measured 
(d, top) and computed (d, bottom) rotational temperature map. Image reprinted from Ref. 151 with 
permission of the authors and the publisher. 







Ref. 151, shows measured rotational temperature compared with computation, shown in the 
lower half of Fig. 2.10(d). Stated root-mean-square (RMS) uncertainties in the velocity 
measurement were -5% with high signal-to-noise. 151 The stated RMS uncertainties in rotational 
temperature ranged from 9% to 35% prior to the Mach disk. 151 

NCF-NO flow tagging velocimetry has also been used to study boundary layer transition on a 
flat plate in a Mach 10 wind tunnel. 153,154 Pure NO 2 was seeded into the boundary layer through a 
spanwise slot located downstream of the sharp leading edge. The angle of attack of the flat plate 
was 20 degrees, reducing the edge Mach number to about 4.2. Parallel focused beams of 355 nm 
light oriented normal to the surface of the flat plate (and in a plane parallel with the streamwise 
direction of the flow) dissociated the NO 2 , creating NO. After a 40 ns delay, the NO was probed 
by a 226 nm laser, at which time the camera acquired an image. One microsecond later, a second 
226 nm laser probed the NO and the second image was acquired with the same camera. The two 
images were compared to compute velocity profdes. Figure 2.11 shows the resulting velocity 
profdes for two cases, one with no trip (i.e. tripping element, a protuberance designed to trip the 
flow from laminar to turbulent) and one with a 1-inm tall cylindrical trip, where the 
measurements were made downstream of, and on the centerline of, the trip. (The laminar 
boundary layer thickness was also approximately 1 mm thick.) The figure shows mean profiles 
(black, top two charts) as well as profiles of the fluctuating component of velocity, u’ (red, 
bottom two charts). For the case of no trip, the velocity profiles compare well with a 
compressible Blasius solution (shown in light grey). When the trip is present, the measured mean 
profiles depart from the laminar solution, showing a profile that is more full than the laminar 
profile near the plate surface with a pronounced velocity deficit near the edge of the profile. The 
fluctuating streamwise velocity increases by a factor of three, up to 250 m/s, between the two 
cases, with the fluctuations highest in the wake of the trip. The single-shot measurement 
precision was 15-25 m/s, which was 1-2% of the maximum velocity in the boundary layer. The 
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accuracy was estimated to be 5-15 m/s. 

2.10 Advantages and Limitations of Laser-induced Fluorescence 

The time scales associated with inducing fluorescence via laser excitation are typically a few 
hundred nanoseconds, which is much shorter than hypersonic flow time scales, therefore 
providing sufficient temporal resolution for high-speed transition-to-turbulence measurements. 
The availability of pulsed lasers capable of picosecond and femtosecond excitation allow for 
fluorescence measurements with time scales much less than those associated with collision and 
reaction time scales. Recently developed kHz- and MHz -rate pulsed laser systems have allowed 
image sequences consisting of tens to thousands 155 of images to be acquired, providing time- 
resolved information pertaining to high-speed fluid dynamic behavior. Both fluorescence-based 
velocimetry 150 and visualization 156 ’ 157 experiments in hypersonic flow fields have been 
perfonned with these types of laser systems. The spatial resolution of a LIF technique is also 
sufficient for many applications, and higher than other techniques such as Raman or CARS, with 
laser sheet thicknesses typically in the range of 0.1 to 1 mm and magnifications of tens of pixels 


per mm, depending on the experimental setup. Fluorescence-based measurements are more 
sensitive than other techniques (for example Raman spectroscopy) with sensitivity on the order 
of parts-per-million or better . 42 
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Figure 2.11: Streamwise velocity profiles on a flat plate in a Mach 10 facility for the case of no trip (a) 
and a k = 1 mm tall, 4 mm diameter cylinder trip (b) and the fluctuating streamwise velocity for no trip 
(c) and a k = 1 mm tall, 4 mm diameter cylinder trip (d). Images reprinted from Ref. 154 with permission 
of the authors. 


Additionally, a wide range of species including intermediate combustion species can be 
probed using fluorescence techniques. Reference 43 provides an extensive listing of many 
species that have been detected using LIF and other methods. Another advantage of LIF is that it 
is readily extended to planar or volumetric measurement (see Refs. 158 and 159). 

Several factors complicate acquisition and interpretation of LIF signals, complicating 
quantitative measurements. Quenching of the fluorescence, which prevents easy quantification of 
signal intensities, was discussed extensively above. Absorption of laser energy as the laser light 
passes through the flowfield can limit the effectiveness of fluorescence-based measurements, as 
the energy will decrease in an exponential manner over a given spatial path length according to 
the Beer-Lambert relations in Eqs. 2.5 and 2.6. This makes quantitative measurements difficult, 
as the laser energy at a particular location may not be easily determined. Absorption can be 
significant when the concentration of the absorbing ground state population is high, the transition 
cross-section is relatively large, the Einstein B coefficient for stimulated absorption is relatively 
large, and/or the path length through which the laser radiation passes is relatively long. To avoid 
strong absorption, a transition may be selected for which the population is small based on 
analysis of the Boltzmann fraction, as was done for PLIF visualization measurements in a 
hypersonic shock tunnel described in Ref. 160 . Absorption can also limit measurement 
capabilities when fluorescence from the probed volume is re-absorbed by the species of interest. 
This effect, known as radiative trapping, occurs when fluorescence emission at frequencies 
readily absorbed by highly populated states must pass through gas containing these potential 
absorbers before reaching the imaging system. Measurements in a non-uniform or turbulent 
mixture are especially susceptible to errors associated with absorption and radiative trapping 
effects, as the absorption coefficient is a spatially- and temporally-varying property. A discussion 
of these issues, and some methods used to circumvent them, is provided in Ref. 42. 

Another disadvantage of LIF is that it usually probes only a single species at a time, 
compared to Raman or CARS, which can interrogate many species simultaneously. However, the 
fact that PLIF can measure spatial distributions of species can compensate for the single-species 
capability in some applications. 

Consideration must also be given to the optical access of the test facility. Since much of the 
work described in this discussion requires fluorescence excitation using laser frequencies in the 
UV portion of the electromagnetic spectrum, the window material used in hypersonic facilities 
must be capable of transmitting such frequencies with minimal absorption. Also, typically two or 
three windows are required for LIF or PLIF applications. The laser (beam or sheet) is typically 
brought in from one window and observed through another window at right angles to the first 
window. A third window can allow the laser to leave the test section, reducing scattered light, 
and allowing the quantification of absorption of the laser beam/sheet in some applications. These 
windows are typically relatively large compared to those required for CARS or diode laser 
absorption measurements. When short-pulse LIF experiments are performed, the inverse 
relationship between the pulse temporal width and spectral width (which for a Gaussian pulse is 


r ia Se r = 0.44/Av) 45 can result in distortion of the temporal characteristics of the pulse. This is 
because the refractive properties of the optical windows of the test facility affect the speed with 
which the frequency components of the short pulse pass through the window material. Reference 
45 discusses issues related to short pulse excitation and interaction with optical components. 

3. Rayleigh and Raman Scattering 

3.1 Introduction 


When a light beam passes through a gaseous medium, it can interact with the gas molecules 
or particles in the gas, thereby scattering light away from the path of the incident beam. Elastic 
scattering occurs if no energy is gained or lost to the medium. If energy is either absorbed or lost 
by the medium, the scattering is inelastic. Light scatter from particles that have a diameter, d, on 
the order of or larger than the light wavelength, A, is termed Mie scattering. Scattering for which 
d« A is known as spontaneous Rayleigh scattering if elastic and spontaneous Raman scattering 
if inelastic. These processes are shown schematically in Figs. 3.1 and 3.2. Represented on the 
energy level diagram in Fig. 3.1, a photon 
excites a molecule from an originating 
state, 1, to a ‘virtual’ state, 2, from which 
the scattered photon is emitted. The 
virtual state, represented by a dashed line, 
is not an actual resonant state of the 
molecule. Instead, it indicates a non- 
resonant, short-lived state in which the 
electron distribution of the molecule is distorted. This virtual state immediately relaxes to the 
originating state (in the case of Rayleigh scattering) or another state (in the case of Raman 
scattering). Relaxation to a higher lying (e.g. vibrational) state than the originating state is 
tenned Stokes Raman scattering. In this case, the molecule absorbs a quantum of energy through 
this process. In anti-Stokes Raman scattering the molecule imparts a quantum of energy to the 
scattered photon so that the scattered photon has higher energy than the incident photon. In this 
case, the originating state must not have been a ground state. 



Figure 3.1. Energy level diagram indicating incident 
radiation, Rayleigh scattering and Raman scattering. 
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Figure 3.2 shows notional Raman/Rayleigh 
spectra, not drawn to scale. Rayleigh scattering 
is shown at the laser’s wavelength. Discrete pure 
rotational Raman lines, associated with rotational 
quanta imparted to or subtracted from the 
incident laser frequency, are shown on opposite 
sides of the Rayleigh peak. Vibrational Raman 


Figure 3.2. Notional Raman/Rayleigh spectra. 

bands are located further away, spectrally shifted towards the red (Stokes) and the blue (anti- 
Stokes). The vibrational Raman bands show rotational fine structure. Raman scattering is much 


weaker than Rayleigh scattering. Raman scattering is typically three orders of magnitude smaller 
than Rayleigh scattering for most gases of interest to supersonic and hypersonic flows. 42 

3.2 Theory of spontaneous Rayleigh and Raman scattering 

Comprehensive reviews of the theory and application of Rayleigh and Raman scattering have 
been given by others. 42,43,161,162 Herein we provide a brief introduction and overview. Later 
sections will show how this theory is applied to measure thermodynamic properties. An 
electromagnetic wave incident upon a molecule perturbs the molecule’s electron cloud, making it 
oscillate at the same frequency as the incident wave. These oscillations cause a periodic charge 
separation within the molecule, known as an induced dipole moment. Oscillating dipole 
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moments act like antennas, emitting radiation. As described in more detail by Baldwin, if 
these oscillations are in-phase, the emission adds constructively, producing a coherent beam. If 
out-of-phase (for example, in the direction orthogonal to the beam), the emitted light interferes 
destructively and the radiation cancels. For a monochromatic plane wave passing through a gas, 
the constructive interference occurs only in the forward direction. The resulting emitted coherent 
light is perfectly in-phase with and indistinguishable from the incident wave. In a gas composed 
of a finite number of molecules, the destructive interference at other angles is not fully complete 
because of statistical variations in the number of particles located in different wavelength-sized 
volumes of the gas.163' 164 That is, there are not exactly the same number of particles in each 
wavelength sized volume, which would be required to cancel out the radiation perfectly. This 
statistical variation in the number density then leads to Rayleigh and Raman scattering. Rayleigh 
and Raman increase significantly in intensity at shorter wavelengths: both scale approximately as 
1/k 4 . 42 

The radiant intensity, 1°, which is the scattered power per unit solid angle, is proportional to 
the square of the induced dipole moment. The induced dipole moment, p, is given by: 42 

p — s 0 a E (3.1) 

where s 0 is the pennittivity of free space, a is the molecular polarizability and E is the incident 
electric field given by E — E 0 cos (<n 0 t), where E 0 is the amplitude of the electric field, u> 0 is the 
frequency of the laser light and t is time. 

The polarizability of a molecule depends on its internal structure and varies with time during 
vibrational oscillations at the natural frequency of the molecule, co v , vibrating in the direction of 
its normal spatial coordinate, Q. The polarizability can be approximated with a Taylor series 
expansion: 

a = a ° + © 0 dQ 


(3.2) 


where the small physical displacement, dQ, of the atoms about their equilibrium positions 
(denoted by the subscript 0) during vibrations is: 



(3-3) 


(3-4) 


(3-5) 


The first term on the right side of Eq. (3.5) describes Rayleigh scattering at a frequency 
corresponding to the incident laser’s wavelength. The second term indicates Raman scattering 
which is shifted from the Rayleigh scattering by ±o> v resulting in upshifted (anti-Stokes) and 
downshifted (Stokes) Raman scattering. While significant theory has been developed to describe 
the physics of Raman and Rayleigh scattering, the strengths of the scattering for different gases 
is generally measured and reported as a temperature-independent differential cross section: 42 


and vary with the laser wavelength but are independent of pressure and temperature. 

Rayleigh scattering from different molecules cannot be distinguished spectrally, so it is not 
usually used to detect individual species. Under circumstances where the composition is fixed or 


measure the gas density, p. As discussed further below, the Rayleigh scattering cross section 
varies from molecule to molecule, with some hydrocarbon species having cross sections more 
than an order of magnitude larger than N?. The cross section for a mixture of gases is equal to 
the mole-fraction-weighted average of the individual cross sections. If the composition can be 


the gas mixture can be determined, allowing the density to be determined from the measured 
Rayleigh scattering. Under conditions where the pressure and composition are kn own or can be 
determined and where the perfect gas law applies, the gas temperature can be inferred from the 



(3.6) 


which can be rearranged as: 



(3.7) 


where the subscript zz refers to a polarization in the z (vertical) direction caused by an incident 
electric field oriented in the z direction, N is the number density of the gas and / is the laser 
irradiance. 42 The differential cross sections vary by process (Rayleigh vs. Raman), by molecule, 


known, or limited to vary under controlled conditions, 165 ’ 166 Rayleigh scattering can be used to 
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estimated or measured, for example by Raman scattering, then the effective cross section for 


42 

measured density. Such measurements can suffer from interference from Mie scattering from 
particles and laser scattered light. To avoid Mie interferences, the gases can be filtered to 
remove particles. Careful experimental design and blackening of surfaces can minimize 
scattered laser light. 


Shift, V 


Spectral analysis of Rayleigh 
scattering yields additional parameters that 
can be measured. Figure 3.3 shows a 
schematic of the dispersed Rayleigh 
scattered light (solid line), separated from 
the incident laser light (dashed line). The 
spectral shift of the Rayleigh from the 
laser is caused by the Doppler shift of the 
gas relative to the incident light beam and 
depends on the detection angle. The width 
of the Rayleigh spectrum depends, in part, 
on the gas temperature. The amplitude of 
the scattering, as described above, depends on the gas density (and composition). Thus, in 
principle, by spectrally dispersing the Rayleigh scattered light, the temperature, velocity and 
density of a gas can be measured simultaneously. These effects are detailed below with 
examples given. 



Figure 3.3. Notional sketch of spectrally dispersed 
Rayleigh scattered light (solid line) which has been 
Doppler shifted from the incident laser light (dashed). 


Revisiting Eq. (3.5), notice that for Raman scattering to occur, ( must be nonzero. This 

happens when a molecule exhibits a change in its polarizability with vibrational displacement. 
For example, as the atoms in the N 2 molecule vibrate, they become less polarizable as the atoms 
approach each other (charges exhibiting more force because they are closer to each other) and 
more polarizable when further apart (charges more weakly interacting). This results in a nonzero 


( — ) . Such vibrations are Raman active. On the other hand, some vibrational modes do not 
\dQJ 0 

change the molecule’s polarizability, for example the asymmetric stretch of CO 2 has the 


following vibration pattern: 4r C 0~>. Such modes have — 0 and are termed 

Raman inactive. Further quantum mechanical selection rules and molecular structure 
considerations (atomic mass, bond length, moment of inertia, etc.) determine the shape and 
structure of the Raman spectra. 42 Since the frequency of Raman spectra depend strongly on the 
individual species’ molecular structure (each vibrational resonance occurring at a different 
energy), Raman spectra from different molecules appear spectrally separated when they are 
dispersed, typically using a grating-based spectrometer equipped with a camera to acquire the 
spectra. Because Raman scattering is species-specific, it can be used to measure individual 
species concentrations, where in Eq. (3.6), N is the species being detected and the differential 
cross section is that of the individual species. 


3.3 Advantages and Disadvantaged of Rayleigh and Raman Scattering 

Rayleigh and Raman scattering have many inherent advantages compared to other 
measurement techniques. Both use just a single laser and the laser need not be resonant with any 
particular molecular resonance as in absorption or laser induced fluorescence (LIF). This allows 
high-powered, fixed frequency lasers to be used for Rayleigh or Raman. Since they are linear 
laser techniques involving only a single excitation beam, Raman and Rayleigh are relatively easy 
to set up and understand, and the data are, in general, easier to analyze than nonlinear techniques 
like coherent anti-Stokes Raman spectroscopy (CARS). Neither Raman nor Rayleigh is sensitive 
to collisional quenching, a phenomenon which complicates the interpretation of LIF signals. 
Absolute intensity calibration of both techniques is relatively straight forward and easily 
perfonned. Raman and Rayleigh can be performed simultaneously using the same laser to 
measure species concentrations (Raman), temperature (Raman and/or Rayleigh), density (Raman 
and/or Rayleigh) and velocity (Rayleigh). 

The major disadvantage of Raman and Rayleigh scattering is the low signal intensity. The 
low signal from Raman scattering generally prevents minor species (less than a few percent by 
mole fraction) from being detected using this technique. Because the signals are so low, large 
(low /-number) collection optics are usually used. Ideally, the detection optics need to be placed 
close to the measurement region and, in ducted flows, large windows are required. Hypersonic 
and combusting flows can sometimes be luminous. Natural luminosity, like spontaneous 
Rayleigh and Raman scattering, emits light in all directions. In order to improve the signal-to- 
noise ratio, temporal, spatial, and spectral filtering can be used to collect the desired radiation 
and block unwanted natural luminosity. For example, gated detection of pulsed signals (or lock- 
in detection of continuous signals) can be used. 

3.4 Translational Temperature, Velocity, and Density Measurements 

The gas temperature and velocity can be detennined from the Rayleigh spectrum by 
resolving the Doppler broadening and Doppler shift, respectively. The Rayleigh scattering 
linewidth is typically in the range of 1-6 GHz (0.03-0.17 cm' 1 ) for supersonic and hypersonic 
flow experiments in the range of a few hundred to a few thousand kelvins, depending on the 
temperature and angle of incidence of detection and collection of the light (see Ref. 161 for 
detailed information about the Rayleigh scattering lineshape and angular dependence). A high- 
resolution laser can be used so that the broadening caused by the laser’s lineshape is negligible 
compared to this Doppler broadening. For example, injection seeded, pulsed Nd:YAG lasers 
typically have linewidths of about 0.1 GHz (0.004 cm' 1 ) while continuous sources can have 
linewidths that are orders of magnitude smaller. To spectrally resolve Rayleigh scattering, a 
high-spectral-resolution instrument is required. Two methods are typically used: gas vapor cells 
and Fabry-Perot etalons. 


A gaseous cell filled with iodine vapor can be used in combination with excitation in the 
visible wavelength range, where h has many absorption lines, by placing a low-pressure gas cell 
containing crystalline iodine and E vapor in front of a detector or camera. The diameter of the 
cell should be larger than the diameter of the collection lens. The cell length, pressure and 
temperature of the h cell are chosen to control the L gas concentration and spectral line shape 
and therefore the absorption magnitude and profile. Attenuation of transmitted light by a factor 
of 10 5 can be achieved using high-resolution cw lasers, though it is more difficult to achieve 
attenuations greater than ~10 with off-the-shelf injection-seeded, pulsed Nd:YAG lasers. The 
absorptions can be used to reject spurious scattered laser light, while passing the pressure- and 
Doppler-broadened and Doppler-shifted Rayleigh scattering, allowing background-scatter-free 
density measurements. These sharp absorption features can also be used for temperature and 
velocity measurements by using the steep edge of the absorption spectral profile to provide 
spectral resolution. Since the absorption spectrum of the gas filter is well known, if the laser 
frequency is scanned across the filter, then the Doppler-broadened and Doppler-shifted Rayleigh- 
scattered light will transmit through the filter, being acquired by the camera. Each pixel on the 
camera will then have obtained a convolution of the Rayleigh-scattered light with the absorption 
spectrum. These spectra can then be deconvolved to determine the gas temperature through the 
thermal broadening and gas velocity through the observed Doppler shift. Miles et al. 161 
demonstrated this approach in a Mach 2 pressure-matched jet flow, shown in Fig. 3.4. 
Variations of this technique can provide both time-averaged and single-shot measurements. For 
single-shot measurements, the dynamic range of velocities that can be measured is limited, 
although by increasing the gas cell buffer gas pressure, the dynamic range has been increased by 
Elliot et al,. who made single-shot measurements with -9% uncertainty over a velocity range 
from 200 to 600 m/s. 169 ’ 170 
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Figure 3.4. Temperature (left), pressure (center) and velocity (right) measured in a Mach 2 supersonic 
jet flow using Rayleigh scattering observed through a gas vapor cell using a frequency-scanned, 
injection-seeded Nd:YAG laser. 161 © IOP Publishing. Reproduced with permission of the authors and of 
10P Publishing. All rights reserved. 


By using a Fabry-Perot etalon, single-shot Rayleigh spectra can be obtained, allowing 
instantaneous (in -100 ns) and simultaneous measurement of temperature, density and velocity. 
A Fabry-Perot etalon consists of two planar, reflective surfaces that cause interference, 
dispersing the transmitted light spatially, so that in the focal plane of the etalon, the spectrum of 






the light is separated spatially and can be resolved spectrally. Two different strategies have been 
developed to acquire these spectra. The first uses a CCD camera in the focal plane of the etalon 
to capture the dispersed spectrum. 171 ' 174 The second method uses spatial masks or mirrors to 
direct the Rayleigh-scattered light to single-point detectors such as photo-multiplier tubes 
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(PMTs). ’ Typically, CCD cameras read out more slowly than single-point detectors such as 
PMTs, so using PMTs generally results in higher-speed detection. Measurement rates up to 32 
kHz have been reported using this approach. 176 However, using a CCD camera offers several 
measurement advantages, described below. 

Figure 3.5 shows one such example of CCD-based detection of Rayleigh-scattered light from 
Bivolaru et al. A pulsed Nd:YAG laser is focused into a heated Mach 1.6 jet flow. The 
Rayleigh-scattered light is collected at right angles by a lens system which down-collimates the 
collected light and passes it through a solid etalon. The etalon-processed light is then focused on 
an electron-multiplying charged-coupled device camera (EMCCD) where the interference 
fringes are realized. This image is called an interferogram. An optical fiber (OF) directs some of 
the spectrally narrow laser light into the optical path so that a circular reference fringe will 
appear in the interferogram, as in Fig. 3.6(a). The pair of oval patterns located in the white 
rectangle show Doppler-shifted, Doppler-broadened, Rayleigh-scattered light collected by the 
lenses and dispersed by the etalon. A similar pattern is shown on the left side of the 
interferogram. These four oval patterns originate from four different spatial locations a few mm 
apart in the flow, thus resulting in four simultaneous, spatially-separated measurements. The gas 
velocity, V, and the Doppler shift, A / , are related by: A / = {(k s -k 0 ) ■ V)/X, where X is the 
wavelength of the incident light, k Q is the wave vector of the incident light and k s is the wave 
vector of the collected light. (The wave vector points in the direction of light propagation and has 
a magnitude of 2n/X.) In the vector diagram shown in Fig. 3.5, k ol is the wave vector of the 
incident light and k sl is the wave vector of the collected light. The observed Doppler shift 
measures the velocity component Vj in the direction defined by k sl - k ol , which bisects the angle 
ft between the laser and collection wave vectors. Similarly, a mirror M r reflects this incident 
beam back through the lens, L 2 , and into the measurement volume with wave vector k „2 resulting 
in collected Rayleigh signal having a wave vector k S 2 . This signal is sensitive to the velocity 
component V 2 . The reflected beam was slightly misaligned in the downward direction so that 
the two measurements would be spatially separated on the interferogram, shown in Fig. 3.6(a). 
Furthermore, the geometry of this experiment was constructed so that the Vj would be 
orthogonal to the jet axis, measuring a radial velocity component while V 2 was parallel to the jet 
axis, measuring the axial velocity component. In Fig. 3.6(a), four spatially separated 
measurements of axial velocity are collected in the bottom half of the interferogram while four 
measurements of radial velocity are collected in the top half. Figure 3.6(b) shows the two boxed 
peaks after they have been processed to linearize the interferograms and to bin the data into 
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single spectra. ’ These spectra were fitted with Gaussians to determine the Doppler shift and 
broadening associated with Rayleigh scattering, relative to the reference peaks. Subsequent 


work by this team and others has been able to simultaneously determine the gas velocity, 
temperature and density from similar spectra. 177,178 



Figure 3.5. Two-component interferometric Rayleigh scattering system from Ref. 173. Mirrors are 
denoted by M, dichroic mirrors by DM and lenses by L. PBS is a polarizing beam splitter while NBF is a 
narrowband filter. 1R is an iris and FPC is a focal plane camera used to monitor the beam alignment. In 
this experiment, the dichroic mirrors were required to filter out other laser beams associated with a dual- 
pump CARS measurement that was being performed simultaneously. Reprinted with permission of the 
authors. 

In Ref. 173, Bivolaru et al. reported velocity measurements with a precision of —40 m/s in a 
flow with -1200 m/s, or about 3% of the maximum velocity. The dynamic range of the 
instrument was -3000 m/s. (The dynamic range is mainly determined by the thickness of the 
etalon, which sets the free spectral range, or fringe-to-fringe spectral separation.) Thus, 
expressed as a percentage of the dynamic range of the instrument, the measurement precision is 
-1%. To obtain more precise measurements at lower temperatures, the experiment can be 
designed to use higher spectral resolution. 



Figure 3.6. (a) Fabry-Perot interferogram of Rayleigh-scattered light obtained by laser beams 
from two different directions as well as laser-light, resulting in the circular pattern, (b) the 
linearized Rayleigh-scattered spectrum obtained from the boxed region in (a), showing best fits 
to the reference and Doppler-shifted light. 173 Reprinted with permission of the authors. 



The use of the camera-based approach has some advantages over the PMT-based approach. 
First, it can tolerate (and in fact, benefits from) scattering from stationary surfaces in the flow. 
Such scattering, as long as it is not too large, provides a reference frequency to determine the 
Doppler shift. Second, having a laser frequency reference in each interferogram makes the 
system insensitive to variations in the etalon transmission spectrum or the laser wavelength. In 
typical PMT-based experiments, the etalon must be temperature controlled, vibration isolated (if 
it is an air spaced etalon) and the laser frequency must be carefully controlled. Uncontrolled 
drift in either the etalon or laser would result directly in a systematic error in the PMT-based 
approach while it is automatically corrected in the CCD-based approach. A third advantage of 
the CCD-based approach is that it is somewhat more tolerant of scattering from clusters and 
particles in the flow, lessening the need for gases to be filtered. Scattered light from particles 
typically appears in CCD-based interferograms as spatially distinct, circular artifacts. These can 
sometimes be removed by image processing. However, if the scattering is too large either 
spatially or in intensity, it can corrupt the measurement; even in this case, other spatial locations 
in the flow may yield measurements from the same interferogram. As shown in Fig. 3.6, the 
CCD-based detection allows multiple spatial points to be measured simultaneously, allowing 
measurements at adjacent spatial locations to be correlated. Furthennore, with CCD-based 
detection, it is straightforward to measure multiple velocity components with the same 
instrument, also shown in Fig. 3.6. Finally, CCD-based detection has often been combined with 
pulsed-laser excitation, which has two benefits: it makes the measurement instantaneous, 
occurring in ~10 ns, and the signal intensity is much higher, allowing measurement in much 
lower density flows, such as atmospheric-pressure flames and low-pressure gas flows. 
Instantaneous measurements at 1/5 of atmospheric density or lower are possible with this 
approach. 173 



Figure 3.7. Five-hundred-pulse-average Raman spectrum in a high-pressure CH 4 air 
flame. The excitation laser wavelength was 532 nm. A 532 nm filter blocks Mie, Rayleigh 
and spurious laser scattering as well as some of the low-rotational-quantum-number 
rotational Raman lines. A subframe burst gating (SBG) technique was also used to subtract 
background emission from this spectrum. 179 Figure courtesy of and with permission of J. 
Kojima, NASA Glenn Research Center. 


3.5 Rotational and Vibrational Temperature Measurements 



Envelope (low resolution) 


T*500K 

/ \ 
1 ^ 

1 

/ 

/ 

J. ■ ■ 



Raman scattering is sensitive to individual molecular rotational and vibrational transitions. 
Since the population of these lines depends on (actually defines) the temperature, it is possible to 
measure rotational and vibrational temperatures from spontaneous Raman spectra. Figure 3.7 
shows a typical Raman spectrum obtained in a hydrocarbon-air flame at high pressure. A 
pulse-stretcher was used to extend the duration of the 500-mJ, 532-nm pulse by nearly a factor of 
10, thereby lowering the peak laser power to avoid laser-induced breakdown while maintaining 
high pulse energy to yield a sufficiently high signal-to-noise ratio. This spectrum shows 
many temperature and concentration dependent features. The relative heights and shapes of the 
different bands depend on temperature and 
gas concentration. Modeling, calibration 
and analysis of such Raman spectra can 
yield rotational and vibrational 
temperatures as well as concentrations. 

Several different strategies for temperature 
measurement based on spontaneous 
Raman scattering exist. Rotational 
temperatures can be measured from pure- 
rotational Raman scattering either using 
high-resolution or low-resolution 
detection, as indicated in Fig. 3.8. The 
advantage of low-resolution detection is 
that the spectrum can be acquired 
simultaneously with the same instrument 
used to acquire multi-species spectra like 
that shown in Fig. 3.7. A disadvantage of 
this low-resolution technique is that many 

different species have similar pure rotational Raman shifts, so they overlap in the same spectral 
region close the excitation laser. Such interferences from different species can lead to 
measurement errors. This technique works over a wide temperature range, including at room 
temperature. Alternately, rotational temperatures could be determined by resolving the 
rotational-vibrational Raman scattering, typically of N 2 , " as is often done for CARS. When 
flows are in rotational-vibrational equilibrium, it is more common to measure the temperature 
using vibrational bands, as described below, because they result in higher signal-to-noise ratios 
and consequently, more precise temperature measurements. 
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Figure 3.8. 

rotational Raman spectra 
temperatures and two different spectral resolutions, 
adapted from Ref. 181.0 10P Publishing. Reproduced 
with permission of the authors and of IOP Publishing. 
All rights reserved. 


Vibrational temperatures can be measured from the relative heights of different vibrational 
Raman bands. A commonly used method compares the ratio of the Stokes to anti-Stokes (S-AS) 
vibrational Raman bands of N 2 . ’ N 2 is often used because it is present in high 

concentration in many supersonic, hypersonic and combusting flows, resulting in adequate 


signal-to-noise ratio. Also, N 2 is used because it is relatively well understood and well resolved 
spectrally from other species. Recall from Fig. 3.1 that the anti-Stokes scattering originates from 
an excited vibrational state whereas the Stokes scattering can originate from the ground 
vibrational state. Thus, the integrated intensity of the anti-Stokes N 2 spectrum at 473 nm in Fig. 
3.7 can be compared to the Stokes scattering at 607 nm to determine the temperature. This ratio 
is a monotonic function of temperature. Since this method integrates the signal intensity in the 
two different bands, it improves the signal-to-noise ratio, allowing single-shot determination of 
temperature. Single shot precisions of -20% of the measured temperature have been obtained at 

i sa 

several points along a line using this method. This technique becomes insensitive to 
temperature below about 700 K because the low population in the excited vibrational state causes 
low signal-to-noise ratio. 42 

A second method of measuring the vibrational temperature is to spectrally resolve the 
different vibrational bands of a molecule such as N 2 , though the individual rotational lines need 
not be resolved. The relative heights of the vibrational levels can be plotted on a Boltzmann plot 
or can be fitted spectrally to determine the vibrational temperature. 187 For example, Sharma and 
coworkers used a KrF excimer laser operating at 248 nm to generate spontaneous Raman 
scattering in the Electric Arc Shock Tube (EAST) facility at NASA Ames Research Center. The 
EAST facility, normally operated as a shock tube, was fitted with a two-dimensional converging- 
diverging nozzle so that a shock reflection generated a high pressure (100 atm), high temperature 
(5600 K) reservoir of N 2 gas that expanded through the nozzle. The sudden expansion through 
the nozzle results in vibrational-rotational nonequilibrium, which was monitored with Raman 
scattering. Measurements were performed at different distances downstream to study the 
evolution of the vibrational relaxation. Figure 3.9 shows a sample spectrum generated and 
spectrally fit to detennine the rotational and vibrational temperatures. The vibrational 
temperatures were then plotted versus distance downstream in the flow so that different 
vibrational relaxation models could be tested. Though not explicitly stated in the paper, the 
temperature measurement precision appears to be about 5%. 



Figure 3.9. An experimental Raman spectrum of N 2 , fitted for rotational and vibrational 
temperature (left), and resulting vibrational temperatures measured as a function of distance 
downstream of the nozzle (right), adapted from Ref. 187 with permission of the authors. 


3.6 Species Concentration Measurement 

Rayleigh scattering has been used in some specialized experiments for determining species 
concentration. For example, the Rayleigh scattering cross section for propane is 13.5 times 
larger than for air, allowing mixing to be quantified in propane/air jets. However, Raman 
scattering is much more commonly used to measure species concentrations. Sandia National 
Laboratories have extensively used Raman scattering for major-species concentration 
measurements (along with Rayleigh and LIF measurements for temperature and minor-species 
concentration measurements, respectively ). 91 This system uses a series of four frequency- 
doubled Nd:YAG lasers to produce temporally-stretched 1.8 Joule pulses which are focused to a 
0.22 mm (1/e ) diameter spot size. A separate laser excites LIF of CO. An imaging system 
directs the collected Raman/Rayleigh/LIF scattered light into the optical analysis system shown 
in the bottom-right panel of Fig. 3.10. The system contains two mechanical chopper wheels (one 
“slow” and one “fast”) which are synchronized with the laser to transmit the signals while 
rejecting flow luminosity. It also has a transmission grating and multiple lenses, beam splitters 
and cooled (low noise) CCD cameras. This system images a 6 -mm-long segment of the 
Rayleigh, Raman and LIF probe volume, allowing concentrations and temperature to be obtained 
along a line. The spatial resolution for the Raman measurements was ~10 pixels / mm along the 
6 -mm-long probe volume, though the actual spatial resolution was slightly worse due to optical 
distortions from the flame . 91 

A sample of the resulting concentration measurements are shown in Fig. 3.10. Raman has 
been used to measure No, O 2 , CH 4 , CO 2 , H 2 O and H 2 while CO was measured with two-photon 
LIF excited at 230.1 nm. Temperature was measured from the intensity of the Rayleigh- 
scattered light, using the ideal gas law to convert density to temperature as described above. The 
data are graphed versus temperature to show the correlation of species with temperature and to 
allow comparison with an unstrained flame calculation, shown in blue. The individual (red) 
data points correspond to different single pulses of the laser. Thus, these measurements are 
‘single shot’ having been obtained with flow freezing (-100 ns) time resolution. Data such as 
these have been acquired at many locations in a variety of different flames, allowing the 
development of a large database of turbulent flames. The accuracy of the Raman concentration 
measurements varies from 2 % to 10 % of the measured concentration, while the precision (based 
on one standard deviation) varies from 0.7% to 7.5% depending on the species. The accuracy of 
the temperature measurement, based on Rayleigh scattering signal intensity, was reported to be 
2% with a 1-a precision of 0.75%. 
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Figure 3.10. Scatter plots of mole fraction and equivalence ratio (()>) data at one location in an 
atmospheric pressure, methane-air flame compared with a calculation, and a schematic of the optical 
components associated with the collection and analysis of light from the Sandia Raman/Rayleigh/LIF 
imaging system (bottom right). Adapted from Ref. 91 with permission of the authors and the The 
Combustion Institute. 


4. Coherent Anti-Stokes Raman Spectroscopy 

4.1 Introduction 


Signal 


Comprehensive reviews of the theory 
and application of Coherent anti-Stokes 
Raman spectroscopy, often called coherent 
anti-Stokes Raman scattering, or CARS, in 
a gas have been given by others. 42,43 
Herein we provide a brief introduction and 
overview. CARS is a non-linear optical 
process in which three laser beams interact 
with the gas generating a fourth, laser-like 
signal beam. The energy level diagram of 
this process (each arrow represents a change in state due to photon absorption, an up arrow, or 
emission, a down arrow, and the length of the arrow is proportional to energy change or, 
equivalently, photon frequency) is shown in Fig. 4.1 The pump beam and the lower frequency 



Figure 4.1. Energy level diagram showing the CARS 
process 


Stokes beam interact with the gas, a pump frequency photon is absorbed and a Stokes frequency 
photon is emitted coherently with the Stokes beam, and the gas is excited to a higher energy state 
via a state that in most CARS setups is a virtual state (although it could also be a real state). The 
difference between the frequencies of the upper and lower state, o> pump — (jo stokes , is called the 
Raman shift. A probe beam photon is coherently scattered from this excited state, shifted up in 
frequency by the Raman shift to form part of the signal beam, and the gas molecule returns to its 
original state. Total momentum as well as energy is conserved. Therefore, the momentum of the 
scattered photons equals that of the incident photons, leading to the following equation: 


kpump T k probe k 


Stokes 


+ k 


signal 


(4.1) 


Since the k’s are the wave vectors (length proportional to photon frequency with direction the 
same as that of the beam), this equation allows the direction of the signal beam to be found from 
the directions of the pump, Stokes, and probe beams. 

A physical interpretation of this process is that the interaction of the pump and Stokes beams 
establishes a optical fringe pattern in the gas. If the pump and Stokes frequencies are the same, 
as with a related technique called degenerate four-wave mixing (DFWM), this fringe pattern is 
stationary. In CARS, the frequencies are different and the fringe pattern moves across the 
interaction region, modulating the intensity at any particular point at the beat frequency, i.e., the 
Raman shift frequency. This moving fringe pattern excites a polarization response in the gas 
which acts as a moving grating from which a probe beam photon is scattered, in a manner similar 
to Bragg scattering, to form the signal beam photon. Since the grating is moving, the frequency 
of the signal photon is shifted relative to the probe photon by the Raman shift frequency. (An 
analogous effect occurs in an acousto-optical modulator where acoustical waves in a solid 
material, such as glass, form a moving grating of varying index of refraction, and incident laser 
light is coherently or Bragg scattered from this volumetric “grating” with its frequency shifted by 
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the acoustical frequency. ) 
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Figure 4.2. Planar BOXCARS beam geometry 


If the pump, probe, and Stokes lasers 
are all single frequencies, then the signal is 
also at a single frequency, as indicated in 
Fig. 4.1. If the pump and/or Stokes lasers 
are broadband lasers, while the probe is 
single frequency, then the signal is 

broadband also, and contains a spectrum that reflects the variation of CARS susceptibility of the 
molecules in the probe volume as a function of Raman shift. In many broadband CARS setups, 
the probe and pump frequencies are the same (derived from the same laser source). In the dual- 
pump CARS technique, 42 ,19 ° pump and probe frequencies are different and these two laser 
beams have interchangeable roles (each beam performs the role of pump in one CARS process 
and probe in a second CARS process). The signal is thus the coherent superposition of the signal 
from the two processes, generated over two different ranges of Raman shift. When the desired 


Raman shift is small, as in rotational CARS (i.e., when the transitions probed are pure rotational 
transitions, which differ by relatively small quanta of energy and are thus relatively closely 
spaced spectrally), the pump and Stokes beams can be derived from the same broadband laser 
while the probe beam is single frequency (dual-broadband CARS). 42 

Since CARS is a non-linear process (see Section 1.2) it requires high irradiance levels. 
Consequently, CARS signal in a gas is typically generated by focusing pump, Stokes, and probe 
beams, and overlapping them at their common focus. If the beams are initially separate and 
crossed at the common focus, the arrangement is called BOXCARS; if they lie in a plane (Fig. 
4.2), it is called planar BOXCARS. In this arrangement the measurement volume is the small 
region of overlap of all three beams. 

4.2 CARS Theory 42 ’ 43 189 

CARS, like Raman and Rayleigh scattering, arises because of time-varying polarization 
induced in the gas in the presence of electromagnetic radiation. CARS, specifically, arises due to 
the third order susceptibility, Xcars > for which the induced polarization is described as follows: 

^(S^signal) — £ oX CARsE (.^ pump) E (.^Stokes) E (.^probe) (4.2) 


The E ' s are the complex electrical field amplitudes. This equation can be substituted into the 
wave equation relating the electrical field to the induced polarization and solved by integration 
along the direction of the signal beam. The CARS irradiance is thus: 
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L is the length of the measurement volume, i.e., the length of the region along which the pump, 
Stokes, and probe beams all overlap. Rigorous calculations of the CARS susceptibility require 
quantum mechanical treatments, but classical derivations are simpler to understand. The gas is 
modeled as a simple harmonic oscillator in which the variation of the normal coordinate, Q 
(intermolecular spacing in the case of vibrational states), with time is described by a second 
order linear ordinary differential equation, with a time-dependent forcing function proportional 
to the average over an optical cycle of the square of the electrical field. (Thus, the forcing term 
oscillates at the Raman shift frequency.) The constant of proportionality in the forcing term is 


proportional to , and consequently to the Raman scattering cross section, where a is the 

optical polarizability of the molecule. The CARS susceptibility is found by solving for the time 
dependence of the normal coordinate, combining with Eqs. (4.1) and (4.2) for the polarization 
(from the chapter on Rayleigh and Raman scattering), and comparing with Eq. (4.2): 
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The susceptibility is proportional to the molecule number density, N, and the differential Raman 
cross section, and reaches a peak where the detuning, A co = co v — ( co purnp — (^stokes) ■> 

approaches zero. The damping coefficient, T, determines the line width. A quantum-mechanical 
treatment yields the following equation: 

+ Xnr where Kj oc NAj Q . (4.5) 

The index j refers to a particular transition and A ; is the fractional difference in the population of 
the two states of the gas between which a transition is taking place (the two real levels in the 
energy level diagram of Fig. 4.1); Xnr is the non-resonant susceptibility, a nearly constant term 
for the CARS process in which all the states in Fig. 4.1, other than the ground state, are virtual 
states. The (real part of the) line shape in Eqn. (5) is a Lorentz function with full-width at half 
maximum equal to r ; - . This line shape is the same as the (complex conjugate of the) line shape 
for the classical solution when the detuning is small relative to the vibrational frequency. The 
line width depends inversely upon the lifetimes of the two states involved in the transition, which 
depend upon the rate of collisions between molecules (collisional broadening). Light emitted by 
molecules in motion is observed at a slightly different frequency by a stationary observer by an 
amount proportional to the velocity component towards the observer. This effect, when averaged 
over many molecules travelling in different directions, results in Doppler broadening. Since the 
distribution of molecular velocity components about the bulk mean in a gas in equilibrium is 
Gaussian, the line shape becomes a Voigt profde, a convolution between Gaussian and 
Lorentzian functions. At pressures much above atmospheric, additional line narrowing effects 
occur. 

Since the CARS signal is proportional to \Xcars \ 2 > a CARS spectrum reflects the 
populations of the molecular energy states involved in the transitions. Since these populations, in 
equilibrium, are related via the Boltzmann equation to the temperature, CARS can measure 
temperature. Integrated CARS signal is also strongly dependent upon number density and, in 
principle, could be used to measure density. This is not typically done because of experimental 
difficulties in maintaining a consistent geometry of the laser beams at the beam intersection. 
Small uncontrolled motions of the beams due to refraction in inhomogeneous gas fields, 
movements of the optical system, variation in the quality of the laser beams, etc., cause changes 
in signal intensity, and calibrations fail. However, where two gas species are resonant in a 
spectrum, the ratio of the population of one species to another may be found from the shape of 
the spectrum. Where only one species present is resonant, but Xnr is known, the fraction of that 
species may be found from the relative amplitude of the resonant signal to the non-resonant 
“background”; however, Xnr depends on the number density of all species present. If all species 
but one are resonant in the spectrum then the composition can be fully determined by reference 
to the non-resonant background. 191 Use of CARS as a diagnostic tool requires comparison of 
experimental and theoretical spectra. Calculation of theoretical spectra is quite complex and 


numerically time-consuming: after computation of theoretical susceptibility, spectra must be 
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convolved with laser line shapes and instrument probe function. 

Nitrogen is a very useful species for 
measurement of temperature since it is 
usually present in fuel-air combustion or 
hypersonic flows, and is readily probed by 
CARS. The band head of the No Q branch is 
located at a Raman shift of 2330 cm' 1 (units 
of inverse wavelength, proportional to 
frequency divided by the speed of light), 
and is readily accessible using available 
lasers. The Q branch occurs as a result of 
transitions between adjacent vibrational 
states (Av=+1, where v is the vibrational 
quantum number) with no change in the 
rotational state (AJ=0, where J is the 
rotational quantum number). Figure 4.3 shows a portion of the N 2 Q-branch spectrum near the 
band head. Transitions are between v=0 and v=land each peak corresponds to a different J value. 
Since rotation of the molecules affects the energy associated with a given vibrational state, the 
energy difference between v=0 and v=l (at given rotational level) is a function of the rotational 
level. Two spectra are shown in this figure - the first is the CARS susceptibility and the second 
is a computed typical broadband CARS spectrum fonned by convolution of typical laser line 
shapes and a typical instrument function with the susceptibility. As may be seen, line shapes 
are typically not fully resolved in experiments. 

Figure 4.4 shows the sensitivity of the signal spectrum to (a) temperature and (b) 
concentration. The calculations are for (a) air and (b) either air (79% N 2 ) or 10% N 2 , 21% O 2 , 
69% H 2 O, with a pressure of 1 atmosphere. The CARS signal strength varies strongly with 
temperature, partly through the effect of temperature on density (N) via the A" dependence of the 
CARS susceptibility. By using the Sandia CARSFT code, the pressure and temperature 
dependence at moderate to high pressures was detennined to be proportional to 194 p 2 T~ 3 5 . 
However, the shape of the spectrum also changes. At low temperature, only one vibrational band 
is present, associated with v=0 to v=l transitions, and is relatively narrow (because fewer J levels 
are occupied). At the higher temperatures, “hot” bands are present also, associated with v=l to 
v=2 and v=2 to v=3 transitions, and more J levels are populated. As the fraction of N 2 is reduced, 
the amplitude of the band structure becomes smaller relative to the non-resonant background, 
and sensitivity to N 2 is lost for concentrations less than a few percent. This lack of sensitivity 
below a few percent is typical and prevents measurement of minor species, although a number of 
CARS techniques, for example resonance CARS (e.g., Ref. 195), have been developed to 
circumvent this limitation. 
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Figure 4.3. CARS N 2 Q-branch susceptibility and 
spectrum near bandhead at 1 atm., 1400 K. 
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Figure 4.4. Effects of (a) temperature and (b) concentration on the CARS N 2 Q-branch spectrum. 


4.3 Ultrafast CARS 

Ultra-fast techniques have been surveyed by Roy et al. 196 The CARS technique as described 
above is a steady-state process involving the mixing of four beams simultaneously present. For 
given laser energy available in the pump, probe, and Stokes laser beams, the energy in the signal 
beam is proportional to t , where t is the time over which the measurement is made. 
Nanosecond (ns= 10' 9 s) CARS is performed with Q-switch pulsed lasers where the pulse energy 
is distributed over 5 ns to 10 ns, at a pulse rate on the order of 10 Hz. There are many 10’s or 
100’s of thousands of Raman frequency cycles during each laser pulse, so that CARS is 
effectively a steady-state process. However, with the advent of femtosecond (fs= 10' 15 s) lasers 
(e.g., titanium-sapphire), very short pulse widths of order 10' 14 seconds are possible at high pulse 
repetition rate (1 kHz or more), and the pulse width is less than the period of a Raman cycle. 
Due to the t scaling, strong signal energy is available with modest pulse energy. 

The theory of fs CARS is different from ns CARS since a steady state problem is replaced 
with an impulsively initiated, time-dependant one. First, a polarization grating is established 
in the gas by the interaction of pump and Stokes beams. This grating evolves in time through 
rotational-vibrational relaxation and then, after some delay (typically up to several hundred 
picoseconds in measurement applications), a signal is generated by scattering the probe pulse off 
the grating. By making measurements at different probe delays the relaxation of the grating, 
which carries with it an imprint of the vibrational and rotational state of the probed molecules, 
may be observed. “Chirping” techniques have been developed where the probe beam is 
broadened both temporally and spectrally in such a way that the frequency of the probe varies 
with time. The temporal relaxation is thus mapped into frequency space, obtained in a single 
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laser pulse, and may be analyzed with an optical spectrometer. 

Picosecond (ps =10' s) CARS employs mode-locked solid state lasers with pulse lengths 
on the order of 10 ps and has characteristics of both ns and fs CARS. As with fs CARS, 
generation of the grating by the pump and Stokes beams, and scattering of the probe to form the 


signal, are typically separate steps. However, the probe pulse is still relatively long compared 
with the Raman cycle period, and the signal contains an optical spectrum similar to that of ns 
CARS. 199 An advantage of ps CARS over ns CARS is that the ps pulses required to obtain 
adequate signal-to-noise ratio CARS spectra are better suited to pass through commercially 
available optical fibers without damaging the fibers. 200 Both fs and ps CARS have the potential 
advantages over ns CARS of much higher data (pulse) rates and the absence, when the probe is 
delayed relative to the pump-Stokes laser pulses, of non-resonant background effects. Another 
advantage of fs CARS is simpler theoretical modeling, since few or no collisions occur in the 
time of the measurement (negligible collisional effects). 196 

4.4 Advantages and Disadvantages of CARS for probing Supersonic and Hypersonic Flow 

In hypersonic propulsion systems (scramjets), flow velocities are supersonic in inlets and 
nozzles, and may be supersonic or a mixture of subsonic and supersonic in the combustors. In 
combustors, pressures are roughly one atmosphere while temperatures are similar to those of low 
speed combustion. Thermal non-equilibrium (not typically present in low speed combustion) as 
well as chemical non-equilibrium can be present due to the short flow through time (on the order 
of milliseconds). CARS has several advantages in this application. The signal comes as a laser 
beam, which means that it can be collected through a relatively small aperture and may be 
separated from non-coherent interferences by spatial filtering. This is particularly useful in an 
engine combustor where there may be emission from the gases and thermal radiation from the 
wall, and where, for structural reasons, it may not be possible to incorporate large windows. 
CARS measurements are spatially and temporally resolved, with ~1.5 mm long and 50 pm 
diameter measurement volume, and a 10 ns time scale. CARS is able to non-intrusively measure 
local temperature and composition. Since the Q-branch spectrum reflects the rotational- 
vibrational state of the molecule, the populations of these states may be directly determined, 
which is useful when thermal equilibrium (and therefore a single “temperature”) does not exist. 

CARS has been used much less frequently in hypersonic freestreams, where temperatures 
and pressures can be very low. Due to the previously-mentioned p 2 T -3,5 scaling of the CARS 
signal, the signal to noise ratio may be low. At low temperatures only the vibrational ground 
state and far fewer rotational states of molecules are populated. Vibrational CARS (e.g., of the 
No Q-branch) then depends upon measuring high resolution spectra of the rotational structure of 
the “cold” band near the band head, which depends upon the pressure-dependent collisional 
effects on line shape, as well as temperature. Measurements of both pressure and temperature, at 
pressures down to about 0.1 atm., have been made in an underexpanded supersonic jet using this 
approach." If the v>0 bands are populated at low rotational temperature (vibrational non- 
equilibrium), this can be easily measured. Alternatively, pure rotational CARS (Av=0, AJ=±2) is 

sometimes used at low temperatures since the lines have greater separation than in the Q-branch, 
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but no information on the vibrational state is obtained. 


CARS has some disadvantages compared to other measurement techniques. It typically 
measures at only a single point, and so does not typically provide simultaneous information at 
multiple spatial locations. Instead, the probe volume is typically scanned around the flow 
allowing time-averaged spatial properties to be measured. The data rate is low in ns systems 
(order 10 Hz) so that in pulsed hypersonic facilities only one single measurement is obtained per 
facility run. Also, its ~1.5 mm long probe volume can be too large in certain applications where 
it is desired to resolve small length scales such as shock waves, shear layers or turbulent eddies. 
CARS requires optical access on two sides of the flow, which may limit application in some 
facilities. The experimental setup for CARS, involving two or three pulsed laser systems along 
with large spectrometer(s), is relatively complicated and time consuming to set up. Furthermore, 
the nonlinear nature of the theory of CARS complicates the interpretation and analysis of the 
resulting spectra. However, with the combination of accurate and precise temperature and multi- 
species measurement capabilities, CARS is commonly used in many supersonic and hypersonic 
flow applications. 

4.5 Temperature and Mole Fraction Measurement 

Extensive dual-pump CARS measurements have been made in a dual-mode scramjet burning 
hydrogen. The CARS lasers, which were located outside the scramjet lab, consisted of an 
injection-seeded Nd:YAG, frequency doubled to 532 nm, an in-house broad band dye laser 
(Stokes laser) centered around 603 nm with FWHM of 10 nm, and a commercial narrow-band 
dye laser centered around 550.5 nm. Beams were relayed to the experiment via a translation 
system that could move the measurement volume and through special slotted windows in the 
scramjet. Beams were focused and crossed in the scramjet in a planar BOXCARS arrangement 
(Fig. 4.2) to fonn the measurement volume. The signal beam was transmitted out of the scramjet, 
recollimated, separated from one of the pump beams, relayed to a 1 m spectrometer, then imaged 
onto a cooled CCD array with 1340x100 pixels. The measured spectra were background 
subtracted and normalized by a CARS spectrum in argon, which has no resonances, to remove 
the spectrum of the broad-band laser. The resulting spectra were fitted to theory, using the 
Sandia CARSFIT code " to generate the theoretical spectra and an in-house fitting code to do 
the fitting, resulting in temperature standard deviations of -3%. 204 The fitted parameters included 
vibrational temperature of No, a single rotational temperature for all resonant species, and mole 
fractions of N 2 , 0 2 , and H 2 . Figure 4.5 shows some typical (averaged) measured spectra and fits 
to theory. Spectrum (a) was in the freestream air of the flow and is shown with a fitted rotational 
temperature of 923 K and vibrational temperature of 1133 K. N 2 and 0 2 Q-branch resonances 
may be seen. The freestream thermal non-equilibrium was a test facility effect. Spectrum (b) is 
in the combustion plume and shows H 2 rotational (S) lines as well as the aforementioned Q- 
branches, with fitted rotational and vibrational temperatures of 1588 K and 1766 K respectively. 
Figure 4.6 shows contour maps of vibrational temperature in the combustor. The wireframe 
represents the corners of the flow path; a single hydrogen injector is seen at the downstream 
surface of the small ramp, on the top surface of the flow path. Flow enters at Mach 2 from the 


test facility nozzle, and is from left to right. The development of the combustion may be seen: 
combustion is initiated on the top of the plume of hydrogen, near the fuel injector, and wraps 
around and engulfs the plume further downstream. 
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Figure 4.5. Typical fits of theory to experimental dual-pump CARS spectra in a dual-mode scramjet: (a) 
freestream, (b) combustion plume. 



Figure 4.6. Contour plots of CARS-measured mean vibrational temperature in a dual-mode scramjet. 

Simultaneous with the temperature measurements, species concentrations were also 
determined from CARS data. The concentration information was derived from the relative 
intensity of the different resonant features in the CARS spectra, comparing to the intensity of the 
nonresonant background, and considering that the mole fractions must sum to one. 191 Figure 4.7 
shows CARS mole fraction measurements obtained at the same conditions as Fig. 4.6. The mole 
fraction is uniform air in the first plane at the left of each of the figures. Evidence of the cold H 2 
fuel jet is seen in the temperature map (Fig. 4.6) and all three mole fraction maps in the second 
plane, which is located just downstream of the fuel injector (Fig. 4.7). As the fuel jet spreads 
spatially and is consumed by combustion, the presence of No in the center of the downstream 
fuel plume shows evidence of fuel-air mixing. However, the O 2 mole fraction does not track 
exactly with the N 2 because it is reacting with the H 2 . Not until the final plane at the right side of 
each figure, when all the H 2 is consumed, does the O 2 penetrate to the center of the duct, in the 
wake of the fuel plume. 
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Figure 4.7. Contour plots of CARS-measured mean mole fractions of N 2 (top), 0 2 (middle) and H 2 
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(bottom) in a dual-mode scramjet at the same conditions as Fig. 4.6. Reprinted with permission of the 
authors. 
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In summary, CARS is a powerful technique for simultaneously and quantitatively measuring 
multiple flow parameters in gas flows with thermal and/or chemical non-equilibrium. If the gas 
is in thermal equilibrium, CARS can be used to measure temperature; if not it can provide 
information on rotational and vibrational (for rotational-vibrational CARS) state populations, 
allowing separate determinations of rotational and vibrational temperature. It also can be used to 
make quantitative measurements of the relative concentration of the species probed. It is 
spatially and temporally resolved, with a short measurement time (10 ns - 10 ps) and a 
measurement volume that is typically order 1 mm long. Nanosecond CARS is limited by 
available lasers to measurement rates of about 10 Hz, but femtosecond CARS data rates are often 
1 kHz or more. Nanosecond CARS is a well-established technique, picosecond and femtosecond 





CARS in gases is quite recent, due to the advent of suitable lasers, and techniques are still under 
development. The primary limitations of CARS include its complexity (which has meant 
relatively few applications of it to large-scale test facilities), the fact that measurements are 
pointwise and not planar or volumetric, and (for nanosecond CARS) the low data rate. Other 
issues include sensitivity to optical misalignment due to vibrations or refractive index variations 
along the beam path, and lack of spectral modeling for some species. 

5. Other Molecular-Based Measurement Techniques 

5.1 Introduction 

Many different surface-based, probe-based and particle -based instruments have been used to 
study transition and turbulence. However, this manuscript focuses only on molecular-based 
measurement techniques. Raman scattering, Rayleigh scattering, laser-induced fluorescence 
(LIF) and coherent anti-Stokes Raman spectroscopy (CARS) are some of the most commonly 
used techniques, and for this reason have been detailed in the preceding sections. This section 
briefly describes additional molecular-based measurement techniques that have been used, or 
that show some promise for use, in studying transition and turbulence in high speed flows. 

5.2 Interferometry 

Interferometry uses interference of light to measure small changes in optical path length. 
Variation in gas flow properties along the optical path cause the index of refraction to vary, 
which generates these path length changes relative to a reference path. The two paths are then 
aligned onto a one- or two-dimensional detector where interference (or a so-called 
interferogram ) is observed. This method is very sensitive because it can measure optical path 
differences that are fractions of a wavelength of the light from which they are generated. 

In fluid mechanical applications, interferometry is usually used for flow visualization or to 
measure the gas density. The gas density is related to the index of refraction through the 
Gladstone-Dale equation: 


n — 1 — Kp (5.1) 

Where n is the index of refraction, p is the density and K is the Gladstone-Dale constant, which 
depends on the gas and the wavelength of light. The Springer Handbook of Experimental Fluid 
Mechanics recently reviewed applications of various types of interferometry to fluid 
mechanics. These include Mach Zehnder interferometry, holographic interferometry, shearing 
interferometry and others. The different interferometric methods can be sensitive to the density 
or the gradient in density. " As a path-averaged measurement technique, most interferometry 
techniques require experimental methods and/or analysis to determine spatially-resolved 
information. If measurements along multiple lines of sight can be obtained, tomographic 
reconstruction can be used to deconvolve spatial information. Or, in axisymmetric flows, an 


Abel inversion can be used. However, in the general case of 3D turbulent and transitional flows, 
these approaches are challenged by the short spatial and temporal scales of the experiments, 
which require many simultaneous high-resolution views from different angles. While time- 
averaged measurements in such turbulent flow fields are possible, it is more difficult to use these 
methods to resolve small temporal or 
spatial scales. 

Laser differential interferometry (LDI) 
is a sensitive method that has been used to 
detect weak perturbations in low-noise 
wind tunnels that have been developed to 
study hypersonic transition to 
turbulence. " This technique uses two 
loosely focused, parallel beams that pass 
through the wind tunnel flow and are 
recombined on a detector on the far side of 
the wind tunnel. A perturbation 
experienced by one of the beams (for 
example, when a laser-generated hot-spot 
passes by) but not the other can be 
observed by the interference of light on a high-speed detector. The system has a flat bandwidth 
response up to 400 kHz, with 0.1 mm spatial resolution (determined by the beam diameter) in the 
streamwise direction. The system can be calibrated, allowing quantitative phase shifts to be 
determined. However, it integrates disturbances along its optical path, so it does not provide 
spatial resolution along the beams. 



Figure 5.1. Schematic of the FLD1 setup in a wind 
tunnel, showing the laser (L), mirror (M), lens (C), 
polarizer (P), Wollaston prism (W), window (B), 
probe volume (A), photodiode (D) and nozzle (N). 208 
Reprinted with permission of the authors. 
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A different version of LDI called focused 
LDI (FLDI) has recently been implemented at 
CalTech 207 ' 208 for studying transitional flows. It 
maintains or even improves the advantages of the 
LDI systems mentioned above in terms of 
sensitivity, dynamic range and fast time 
response, but it has the important advantage of 
being able to spatially localize the measured 
perturbations. Whereas LDI uses two loosely 
focused and widely separated (~7 mm) laser 
beams, FLDI uses stronger focusing and closer 
spacing. Figure 5.1 shows the FLDI 
experimental setup. The Wollaston prisms have 
the effect of deflecting the incoming beam into 
two beams with perpendicular polarizations. A 



Figure 5.2. Calculated laser beam profdes at 
the focus of the FLDI measurement region for 
the two orthogonally parallel beams (one yellow 
and one blue, with green showing the 
overlapped region). 207 Reprinted with 
permission of the authors. 


lens then directs these beams into the flow where they focus. As shown in Fig. 5.2, the two 
beams (one yellow and one blue), focus in two different spatial locations in the flow (because the 
Wollaston prism causes them to travel at slightly different angles before passing through the 
lens). They then diverge, pass through another lens and Wollaston prism, a polarizer, and onto a 
detector where the two beams interfere. 

The advantage of the FLDI technique is that, 
over most of the optical path, the two beams 
have common optical paths so that both beams 
experience the same disturbances. This is shown 
as the green region in Fig. 5.2. Thus, optical 
disturbances in this region (for example, wall 
boundary layers or disturbances elsewhere in the 
flow) cancel out, leaving sensitivity only to 
disturbances in the focus region. The resulting 
interference signal is then mostly sensitive to the difference in density, A p, between the gas 
located at the two foci. Data are usually reported by referencing to the mean local density, p, 
which must be estimated. The two beams focus to a 0.1 mm diameter and are located 0.35 mm 
apart. The reported spatial resolution is 0.7 mm in the streamwise direction, although the system 
is sensitive to perturbations along the ±10 mm interaction length. The system has a 20 MFIz time 
response, which is the highest frequency response of any measurement technique described in 
this manuscript. The system was used to study small disturbances on a cone, with the laser beams 
oriented tangent to and ~0.5 nun above the surface of a cone model, where the boundary layer 
was slightly larger than 1 mm thick. A typical burst originating from a wave packet travelling in 
the boundary layer is shown in Fig. 5.3. The strongest oscillation occurs at 1.17 MHz. The 
reported accuracy in measuring A pip is -30% while the precision is 0.5% in A pi p, sufficient to 
resolve acoustic instabilities in the flow. Limitations of this technique are its somewhat long 
(-20 mm) probe volume and that it is sensitive to the density difference between two closely 
spaced points instead of the easier-to-compute and more physically intuitive density fluctuation 
about the mean. Nonetheless, it is much more sensitive than the other measurement techniques 
reported in this manuscript. 

5.3 Background Oriented Schlieren (BOS) 

A relatively new measurement technique that is sensitive to density (through index of 
refraction gradients) is known as Background Oriented Schlieren (BOS). 209,210 BOS uses a CCD 
camera equipped with a lens to image a flow placed in front of a background consisting of a 
random dot pattern (or other” ’ random pattern). A reference image of the background 
pattern is acquired with the flow absent and then a second image is acquired with the flow 
present. Refractive index gradients cause the dots to move slightly from one image to the other. 
The pair of dot pattern images is processed using particle-image velocimetry (PIV) software to 
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Figure 5.3. FLDI measured instability wave 
packet in the boundary layer of a cone. 208 
Reprinted with permission of the authors. 



recover the path-integrated light deflections (pixel shifts) caused by the refractive index (density) 
gradients. 

To obtain the highest possible spatial resolution, both the dot pattern and the flow 
disturbance should be in crisp focus. This suggests the use of small lens apertures to obtain a 
large depth of field. To obtain maximum sensitivity, the flow disturbance should be located as 
close to the camera lens as possible. Since it is difficult to achieve a long enough depth of 
field to optimize both of these requirements, in practice, the object is usually put closer to the 
background dot pattern, for example from 0.5 to 0.875 of the way from the camera to the 
background pattern. ’ Also affecting spatial resolution are the camera pixel size, the 
resolution and quality of the background pattern used and the interrogation window size used in 
image processing. 

In a two-dimensional flow, the pixel shifts can be used to directly compute the density 
gradient and thus the density field.' In an axisymmetric flow, the density gradient varies along 
the optical path and additional processing of the data is required to make quantitative 
measurements. In such axisymmetric cases, an Abel or Radon inversion can be used to compute 
the density field, for example in the supersonic flow over a cone at zero degrees angle of 
attack. In non-axisymmetric flows, tomography using multiple camera views is required to 
obtain density information." The multiple views can either be obtained using a single camera 
that observes the flow from different angles at different times, resulting in time-averaged 
measurements,' or by using multiple cameras to take many simultaneous images, providing 
time resolution, which is a current research topic. Since BOS measures gradients and is thus a 
relative measurement technique, the density in one point in the flow must be known or estimated 
to obtain the density field. BOS has the advantages of being a low-cost technique and being able 
to provide both qualitative flow visualization and quantitative measurements for large fields of 
view. One disadvantage is that BOS does not produce real-time visualization and measurement 
results since the raw images require processing. Also, as a tomographic technique using a 
limited number of views, it has much lower spatial resolution than laser-sheet based techniques 
such planar laser induced fluorescence, at least in unsteady flow where it is impractical to have a 
large number of views. 

5.4 Tunable Diode Laser Absorption Spectroscopy (TDLAS) 

Another path-averaged measurement technique used to study supersonic and hypersonic 
flows is Tunable Diode Laser Absorption Spectroscopy (TDLAS). ’ TDLAS measurements 
are path-integrated, so experimental methodologies and post-processing are required if it is vital 
to extract spatially resolved information from the data. Though not necessarily providing 
spatially resolved information in all implementations of the technique, TDLAS has several 
advantages compared to the measurements described in prior sections. Relatively inexpensive, 
small, low-power lasers from the telecommunications industry can be used. These lasers can 
scan at very high rates, providing information at tens or hundreds of kHz, ’ and the lasers are 


typically routed to and from the flowfield using optical fibers, simplifying the experimental 
setup. Technology from the telecommunications industry, such as fiber couplers and 
multiplexers, can be used to split the laser light into many channels and/or to combine many 
wavelengths into a single fiber, allowing multiple species to be probed with multiple colors 
along many paths through the flow. The light from the fibers is coupled into lenses, propagated 
across the flow and captured by lens-coupled multi-mode fibers, separated into wavelengths and 
directed onto detectors. Thus, different species and spectral lines can be probed at different 
locations in a ducted flow, for example, allowing some spatially-resolved concentration and 
temperature information to be determined from the line-of-sight measurements. 


Recent advancements in TDLAS include systems that scan a single laser or multiple lasers 
spatially across the flow and from different directions to allow time averaged tomographic 
reconstruction of temperature and species concentrations in turbulent combusting flowfields or 
a grid of fixed optical paths transmitting across the flow, allowing a low-resolution but time- 
resolved tomographic reconstruction of concentration and temperature. ’ Figure 5.4 (from 
Ref. 220) shows temperature and mole fraction maps obtained in a General Electric J85 gas 
turbine engine. Thirty individual laser beams — 15 in the horizontal direction and 15 in the 
vertical direction, making a square mesh of 225 grid points — were used to probe the flow. The 
laser was rapidly scanned over multiple water vapor absorption lines near 1.35 pm with scans 
occurring at a rate of 50 kHz. While the smaller turbulent length (and time) scales of the flow 
are not resolved with this system, the large-scale spatial (and temporal) distribution of the 
temperature and mole fraction are resolved. 
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Figure 5.4. Temperature (left) and water vapor mole fraction (right) measured at the exit of a General 
Electric J85 engine, with 50 kHz diode laser tomography, adapted from Ref. 220. The paper online has 
li nk s to play animated movie clips of these data. Reprinted with permission of the authors and publisher. 





5.5 Degenerate Four- Wave Mixing (DFWM) 


Like CARS, Degenerate Four-Wave Mixing (DFWM) is a nonlinear, four-wave mixing 
measurement technique that involves the crossing of three laser beams to generate a fourth beam 
that is detected. DFWM’s three input laser beams are all of the same color (degenerate) and are 
usually tuned to an absorption in the atom or molecule being probed. The signal beam from 
DFWM measurements is typically analyzed in the spectral domain using scanned or broadband 
lasers to measure a spectrum of signal intensity versus wavelength. DFWM has been used to 
measure gas concentration, temperature and velocity. 42 Compared to LIF, it has an advantage for 

measuring species concentrations because it is relatively insensitive to collisional quenching, 
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owing to the fact that it is an absorption-based technique (not depending on fluorescence). 
Since each beam is resonant, it is also relatively sensitive and can be used for detecting minor 
species. Since the DFWM signal is a laser-like beam, it is suitable for studying luminous flows 
since background radiation can be spatially filtered with an aperture. DFWM is usually 
implemented as a single-point measurement, though it has been extended to line and 2D imaging 
applications. The main drawbacks of DFWM are that it is more difficult to align and 
implement than LIF, it has much lower spatial resolution than LIF (because of its elongated 
probe volume), it is sensitive to absorption of the laser beams, it is sensitive to saturation of the 
molecular transitions, and it requires very high beam quality to obtain 2D images. 42 DFWM 
temperature measurements in a flame have reported accuracies and precisions more than a factor 
of two worse than CARS " ’ while also having a larger probe volume. 

5.6 Laser Induced Thermal Acoustics (LIT A) 

Another four- wave mixing technique, laser induced thermal acoustics (LIT A), uses two 
pump beams to generate acoustic waves that scatter a probe laser into a signal beam that is 
detected." The two pumps and probe can be the same wavelengths, in which case it is also 
DFWM. Or, the probe beam can be a different wavelength than the pumps, which reduces light 
interferences from the incident, high-powered pump beams. Whereas the signal beam from the 
DFWM techniques described in the previous section is usually analyzed in the spectral domain, 
LITA usually uses the time domain (signal beam intensity versus time) to measure gas velocity, 
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speed of sound (which can be converted to temperature if the gas composition is known) ' 
and gas pressure. " In the time domain, the signal intensity is characterized by a damped 
oscillation, which can have beat frequencies containing the velocity and temperature 
in formation. Velocity and temperature accuracies and precisions on the order of 1 m/s and 1 K, 
respectively, have been reported in air flows. These accuracies and precisions are the best 
reported in this manuscript for these parameters. LITA is usually implemented as a single-point, 
single -velocity-component measurement. The main limitations of LITA are its relative 
complexity and that the measurements are obtained over relatively long probe volumes - on the 
order of 1 to 2 cm - which is much longer than most transitional and turbulent length scales 
occurring in high speed flows. Since the technique is a nonlinear technique, the resulting 
measurements are not simple averages over the probe volume. Instead, as with CARS, there is a 


nonlinear bias that skews the averaging towards higher gas density (lower temperature) 
regions - " which complicates interpretation of the measurements. This problem is more severe 
for LITA, which typically has a probe volume ~10x longer than CARS. 

6. Conclusion 

This manuscript has described numerous molecular-based measurement techniques that are 
potentially applicable to the investigation of high-speed transitional and/or turbulent flows. One 
might ask, “Which technique is best?” or “Which one should I use?” The answer to the first 
question is that none is best. The different techniques measure different properties in different 
ways and have relative merits. The answer to the second, more important question depends on 
many factors. It depends on the measurement requirements as well as the past experience of the 
research team. It depends on the available equipment, the budget for new equipment and the 
time frame available to do the measurement. Some techniques, such as CARS, require a year or 
years of training to leam well enough to apply the technique expertly. Dual-pump CARS 
requires home-built equipment (i.e. a broadband dye laser) as well as two commercially available 
lasers, a spectrometer, a low noise CCD camera, and other equipment. It is therefore relatively 
time consuming and expensive to set up. If a researcher or research team has experience with 
similar technology (for example Nd:YAG lasers, dye lasers, CCD cameras), then the work will 
progress faster. Less complicated (easier to set up and to understand) methods might yield 
usable results in a shorter time frame. These practical issues are certainly important, but a major 
consideration in planning an experiment is the measurement requirements. 

A measurement campaign should begin with an interview of the customer for whom the data 
is being acquired. Answers to the following questions need to be obtained: 

• What parameter(s) need to be measured? 

• Must multiple parameters be obtained simultaneously to determine correlations? 

• What spatial resolution is required? 

• Is imaging required or are single-point or line measurements sufficient? 

• What temporal resolution is required? (e.g. time required for a single measurement) 

• Do measurements need to be time resolved? (e.g. a continuous sequence of data) 

• What accuracy is needed? 

• What precision is needed? 

• What quantity of data is required? 

• When is the data needed? Is instant (real-time) data required? 

• Where in the flow are measurements required? (inflow, exit, near walls, etc.) 

• What type of optical access is available? 

• Can (toxic) seed gases be introduced? Will they influence the properties being measured? 

• What is the ordered priority of the above requirements? 


This list of requirements must be matched up with available measurement technologies such 
as those reported in the sections above. There is rarely a perfect solution for this exercise. 
Instead, there is a compromise between the needs of the customer and the capabilities of the 
instrumentation team. A common set of goals should then be established and agreed upon before 
the research actually starts. This exercise is fruitful for the measurement scientist because it 
often results in a need for a new measurement technique that can be developed if existing 
techniques are not suitable. 

In the case of measurements of transitional flows, the precision, time resolution, acquisition 
rate, and spatial resolution requirements are all important factors. Considering only the precision 
requirement, we see from Table 1 . 1 that freestream perturbations in quiet wind tunnels are on the 
order of 0.01% - 0.2%. Even the most sensitive technique reported herein, FLDI, which can 
measure gas density with a precision of 0.5%,' lacks the sensitivity to monitor quiet tunnel 
freestream disturbances (although it is possible that this instrument could be optimized or 
improved to make such measurements). LITA, which has a demonstrated single-shot precision 
of ~1 m/s' could potentially measure these fluctuations. From Table 1.1, it can be seen that 
conventional tunnels have larger fluctuations, typically on the order of a percent. But to measure 
these fluctuations, sub-percent precision is required. FFDI and possibly FITA would appear to 
be appropriate measurement techniques in these applications, but most of the other techniques 
reviewed herein are insufficient as far as resolving fluctuations is concerned. 

From Table 1.2, the fluctuations in the transitional flow downstream of discrete roughness on 
a flat plate are much higher - on the order of tens of percent. Turbulent boundary layer flow 
(also shown in Table 1.2) also shows tens-of-percent fluctuations. To resolve these fluctuations, 
techniques with precisions on the order of a few percent are suitable. Several of the techniques 
described herein are suitable for this application, for example Rayleigh scattering, FIF molecular 
tagging velocimetry, and CARS thermometry, all of which have precisions around a few percent. 
Still other techniques have precisions that are so poor that they cannot even resolve turbulent 
fluctuations. These techniques still may be useful or measuring mean quantities in turbulent 
flowfields. 

Nonetheless, measurement precision is not the only requirement: spatial resolution and other 
factors must be taken into consideration. For example, though FITA may be precise enough to 
measure freestream perturbations in conventional wind tunnels, it is possible that the 1 cm long 
probe volume may be too large to observe small-scale flow structures. Also, FITA uses high- 
powered lasers and has usually been implemented at a rate of ~10 Hz, which is far too slow to 
provide time-sequence data that may be required in some applications. Turbulent flows, by 
definition, have a wide range of length scales. Many of the measurement techniques described 
herein may be able to resolve the larger length scales, but great care must be taken in order to 
resolve the smallest structures. Even if only larger turbulence length scale data is obtained, the 
data can be of value to the high-speed turbulence research community, especially since very little 


experimental data (both mean and, particularly, fluctuating) have been obtained in the hypersonic 
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regime. 

Clearly, more sensitive (higher precision) measurement techniques are needed to study 
transition and turbulence in high speed flows. It is also desired to make the measurements 
simpler and easier to understand, to set up, and to use. Higher repetition rate and higher 
accuracy methods are also needed. Since no single technique can measure every parameter of 
interest with uncertainties that meet all requirements, combining existing instruments together is 
a possible method for meeting more requirements. For example, CARS/Rayleigh" has been 
used for temperature, major species concentrations and velocity in unseeded supersonic 
combusting flows, Raman/Rayleigh/LIF 91 for major and minor species concentrations and 
temperature in turbulent combusting flows, and CARS/LDV for temperature and velocity in 
particle-seeded combustion flows. However, performing multiple measurements simultaneously 
brings additional complexity, and compromises often must be made between the different 
technologies. An alternate approach is to extend the capabilities of an existing technique to 
measure more parameters. For example, Tedder et al. used WIDECARS to extend the dual- 
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pump CARS technique to measure additional species." 

With continued research, molecular-based instrumentation for studying high-speed 
transitional and turbulent flows will continue to improve in its specific capabilities (e.g., 
precision, accuracy, spatial and temporal resolution), as well as in ease of implementation. Close 
collaboration with customers will help define research goals. Ultimately, the use of molecular- 
based diagnostic techniques will provide data to help validate computational codes as well as to 
help understand the complex fluid physics of high-speed unsteady flows. 
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